5*  8  12 


4Kuftltv  (tLASSiFiCATlON  OF  THi$  PAGE 

|  REPORT  DOCUMENTATION  PAGE 

la  REPORT  SECURITY  CLASSIFICATION 

1b  RESTRICTIVE  MARKINGS 

UNCLASSIFIED 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION /AVAILA8ILITY  OF  REPORT 

Approved  for  public  release; 

2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

Distribution  is  unlimited 

4  PERFORMING  ORGANIZATION  REPORT  NuMBER(S) 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6*  NAME  Of  PERFORMING  ORGANIZATION 

Naval  Postgraduate  Schoo 


6c  ADDRESS  (City.  State,  end  ZIP  Code) 

Monterey,  California  93943-5000 


7a  NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


7b  ADDRESS  (City.  State,  end  HP  Codt) 

Monterey,  California  93943-5000 


8a  NAME  OF  FUNDING  -SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c  ADDRESS  (City,  State  and  HIP  Code ) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO 


11  title  (include  Security  Classification)  C0MPUTER  SIMULATION  MODEL  FOR  STUDYING  AIRCRAFT  TAKE¬ 
OFF  SCHEDULES  AT  A  TRAINING  AIR  FORCE  BASE 


iz  PERSONAL  author(S) 


Macropoulos,  Dimitris  G. 


13a  TYPE  OF  REPORT 

Master's  Thesis 


■3D  T  me  COVERED 
PROM  TO 


14  DATE  OF  REPORT  (Year,  Month,  Day)  IS  PAGE  COUNT 

1988  March  96 


i6  supplementary  notation  ine  views  expressed  m  this  thesis  are  those  of  the  author 
and  do  not  reflect  the  official  policy  or  position  of  the  Department  of 
Defense  or  the  u.S.  Government. 


17  COSATi  CODES _ I  '8  SUBJECT  TERMS  (Continue  on  revtne  if  necessery  end  identify  by  block  number) 

F'ELD  |  GROUP  |  SUB  GROUP  |  m  —  . 

Simulation,  Air  Traffic 


19  ABSTRACT  (Continue  on  reverse  if  necessary  end  identify  by  block  number) 

A  computer  simulation  model  for  studying  take^pff  schedules  at 
Kalamata  Air  Force  Base  in  Greece.  Six  aircraft  take'off  schedules  were 
examined  and  a  comparison  of  results  was  based  upon  performance  and 
® f ^ iciency/saf ety .  The  overall  simulation  model  can  be  easily  modified 
to  examine  other  aircraft  takeoff  schedules,  fc ,  .V 


r'  i-'  ,  C'Ka  -fl 0-\4if 


-I  X'S] 


20  DISTRIBUTION 'AVAILABILITY  OF  A8STRACT  1 21  ABSTRACT  SECURITY  CLASSIFICATION 

CXUNCLASSIFIEO/UNUMITED  □  SAME  AS  RPT  □  DTlC  USERS  I  UNCLASSIFIED 


ZZb  TELEPHONE  (Include  Are*  Code)  tic  OFFICE  SYMBOL 

(408)  646-2413  55As 


DO  FORM  1473, 84  MAR  83  APR  Ration  may  be  used  until  exhausted  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

AM  other  editions  are  obsolete  .  u.».  _ _ „  *,,,«*  e„,„ 


1 


Approved  for  public  release;  distribution  is  unlimited. 


Computer  Simulation  Model  for  Studying  Aircraft  Take-off 
Schedules  at  a  Training  Air  Force  Base 

by 

Dimitris  G.  Macropoulos 
Captain,  Hellenic  Air  Force 
B.S.,  Hellenic  Air  Force  Academy,  1978 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 


Author : 


NAVAL  POSTGRADUATE  SCHOOL 
March  1988 


Dimitris  G.  Macropoulos 


Approved  by: 


A.  F.  Andrus,  Thesis  Advisor 


Peter  Purdue,  Chairman 
Department  o£^Operatl<^ns  Research 


JamesT  1 M.  Fremgea 

Acting'' Dean  of  ^Information  and  Policy  Sciences 


ABSTRACT 


This  thesis  presents  a  computer  simulation  model  for 
studying  take-off  schedules  at  Kalamata  Air  Force  Base  in 
Greece.  Six  aircraft  take-off  schedules  were  examined  and 
a  comparison  of  results  was  based  upon  factors  of 
performance  and  efficiency/safety.  The  overall  simulation 
model  can  be  easily  modified  to  examine  other  aircraft  take¬ 
off  schedules. 
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I .  INTRODUCTION 


THE  PROBLEM 

Air  Traffic  Control  has  become  a  safety  issue  of  great 
importance  during  the  last  decade  because  of  the  many 
near-miss  or  tragic  accidents  that  have  occurred  at  military 
and  civilian  airports  worldwide.  The  main  reasons  for  these 
accidents  are: 

*  Air  Traffic  Control  system  failure, 

*  Air  Traffic  Control  erroneous  procedures, 

*  Pilot  error, 

*  Weather  conditions, 

*  Increased  air  traffic. 

This  issue  has  even  greater  significance  at  military  air 
training  bases  because  of  the  very  high  ratio  of  student 
pilots  to  experienced  pilots  that  are  using  the  air  space 
and  because  of  the  large  volume  of  aircraft  activity  in  the 
air  in  specific  areas.  It  is  standard  military  training 
practice  for  a  large  number  of  training  aircraft  to  be 
assigned  to  the  same  radio  channel  for  aircraft-to-aircraft 
and  for  aircraft-to-air  traffic  controller  communications 
and  to  have  a  large  number  of  training  aircraft  following 
approximately  the  same  air  pattern  with  the  same  air  speed 
and  altitude.  These  conditions  can  increase  the  probability 
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of  breakdowns  in  synchronization  and  communication  between 
aircraft  and  between  the  aircraft  and  the  air  traffic 
controllers  and  can  therefore  increase  the  probability  of 
accidents  occurring. 

Kalamata  Air  Force  Base  is  an  Air  Force  training  base 
located  in  southern  Greece.  This  Base  utilizes  two  types  of 
aircraft  for  training  purposes,  T-37's  and  T-2's.  The  two 
different  types  of  aircraft  used  at  Kalamata  Air  Force  Base 
ordinarily  have  different  flight  schedules,  flight 
capabilities,  and  training  missions  that  would  allow  them  to 
take-off  and  begin  their  training  flights  at  different 
times.  However  because  of  current  operational  constraints  it 
is  common  for  a  group  of  T-2  and  T-37  aircraft  to  complete 
their  missions  at  approximately  the  same  time  so  that  they 
simultaneously  return  to  the  local  Air  Traffic  Pattern.  It 
is  during  the  simultaneous  approach  of  the  returning 
aircraft  that  a  critical  safety  problem  arises  due  to  the 
traffic  congestion  and  pilot  fatigue.  Furthermore,  the 
inefficient  aircraft  schedules  interfere  with  the 
performance  of  the  scheduled  activities.  A  more  efficient 
scheduling  of  aircraft  take-offs  can  remedy  the  safety 
problems  and  also  can  diminish  the  necessity  for  aircraft  to 
wait  in  order  to  get  into  the  mission  areas. 

This  thesis  provides  a  computer  simulation  model 
programmed  in  GPSS  for  the  IBM-PC  for  analyzing  current 
procedures  for  the  efficient  management  and  control  of  the 
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air  traffic  of  Kalamata  Air  Force  Base.  The  thesis  contains 
the  model  logic,  the  GPSS  program,  model  validation  and 
model  results. 
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II.  MODEL  DEVELOPMENT 


A.  DESCRIPTION 

The  environment  that  is  modeled  is  the  operation  of  the 
Kalamata  Air  Force  Training  Base  in  Greece.  This  operation 
is  modeled  as  a  sequential  multiserver  limited  capacity 
queuing  system.  The  server  elements  of  the  model  consist  of 
the  base,  two  runways,  an  air  traffic  pattern  and  eight 
training  mission  areas. 

The  base,  runways  and  air  traffic  pattern  are  illustrated 
in  Figure  2.1.  The  air  traffic  pattern  is  modeled  as  a  set 
of  sequential  servers  consisting  of  two  entrance  points  (EP1 
and  EP2),  the  initial  point  (IP),  the  low  initial  point 
(LIP),  the  break  point  and  the  base  key.  EP1  is  the 
entrance  point  for  the  aircraft  returning  from  a  western 
mission  area  and  EP2  is  the  entrance  point  for  aircraft 
returning  from  an  eastern  mission  area. 

The  mission  areas  are  represented  in  the  model  as  single 
points.  Specific  mission  area  training  activity  is  not 
modeled.  The  eight  mission  areas  are  the  neighboring  areas 
around  the  airport  and  are  illustrated  in  Figure  2.2.  The 
areas  west  of  the  airport  are  numbered  one  through  four  and 
the  areas  east  of  the  airport  are  numbered  five  through 
eight. 
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The  calling  population  consists  of  a  basic  training 
squadron  employing  T-37  type  aircraft  an  advanced  training 
squadron  employing  T-2  type  aircraft  and  occasional 
aircraft  from  other  bases. 

B .  MODEL  LOGIC 

The  general  flow  through  the  model  for  a  typical  aircraft 
is  as  follows: 

A  take-off  time  is  scheduled  for  the  aircraft  from  one 
of  the  two  squadrons  and  the  aircraft  enters  the  base 
activity  queue  at  that  time.  If  there  is  no  runway  or  air 
space  conflict  with  other  aircraft  landing  or  taking-off, 
the  aircraft  takes-off  and  proceeds  to  the  base  departure 
point  where  it  requests  and  is  assigned  to  a  mission  area  to 
carry  out  the  scheduled  training  activities. 

If  all  training  mission  areas  are  occupied  the  aircraft 
is  assigned  to  the  mission  training  area  with  the  smallest 
waiting  line  and  in  the  case  of  ties  the  aircraft  is 
assigned  to  the  mission  area  with  the  earliest  expected 
departure  time  for  the  occupying  aircraft.  All  the  aircraft 
that  wait  for  a  mission  area  to  become  available,  maintain  a 
maximum  altitude  of  7000  feet  at  the  corresponding  areas 
and  perform  training  maneuvers  consistent  with  the  altitude 
safety  requirements.  Once  the  aircraft  enters  the  mission 
training  area  it  stays  in  the  mission  area  for  a  standard 
length  of  time  for  the  aircraft  type  and  performs  the 
scheduled  activities. 


At  the  completion  of  the  mission  area  training 
activities  the  aircraft  returns  to  the  base  and  enters  the 
returning  air  traffic  pattern  and  follows  the  sequence  of 
events  as  described  below  and  illustrated  in  Figure  2.1. 

Aircraft  returning  from  the  western  mission  areas  enter 
the  air  traffic  pattern  at  entrance  point  EP1  and  aircraft 
returning  from  the  eastern  mission  areas  enter  the  air 
traffic  pattern  at  entrance  point  EP2.  The  sequence  of  air 
legs  in  the  base  air  traffic  landing  pattern  and  the 
aircraft  capacity  of  each  air  leg  is  as  follows. 

Air  Leg _ Aircraft  Capacity 


EP1-EP2  2 
EP2-IP  2 
IP-LIP  2 
LIP-Break  Point  3 
Break  Point-Base  Key  3 
Base  Key-Final  2 


If  any  of  the  first  four  air  legs  is  at  full  capacity 
when  an  aircraft  attempts  to  enter,  the  aircraft  must  orbit 
at  the  air  leg  entry  point  and  wait  for  entry.  If  the  Break 
Point-Base  Key  leg  is  at  full  capacity  when  an  aircraft 
attempts  to  enter,  the  aircraft  attempting  entry  returns  to 
the  LIP-Break  Point  air  leg  entry  point.  If  the  Base  Key- 
Final  leg  is  at  full  capacity  when  an  aircraft  attempts  to 
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enter,  the  aircraft  attempting  entry  performs  "go  around", 
that  is,  it  cancels  the  landing  and  attempts  a  "close 
approach"  pattern. 

Aircraft  in  the  "close  approach"  pattern  re-enter  the  air 
traffic  pattern  at  the  Base  Key  point  if  no  other  aircraft 
is  in  the  Break  point-Base  Key  air  leg  and  if  no  other 
aircraft  is  waiting  for  take-off.  Otherwise  they  re-enter 
the  air  traffic  pattern  by  joining  the  waiting  line  at  the 
LIP  entry  point. 

Upon  being  cleared  for  landing  the  aircraft  will  either 
land  and  return  to  the  squadron  or  if  flying  constraints 
allow  the  aircraft  will  perform  a  touch  and  go  landing.  If 
a  touch  and  go  landing  is  made  the  aircraft  will  either 
enter  the  "close  approach"  pattern  or  attempt  to  re-enter 
the  air  traffic  pattern  with  associated  probabilities  of 
0.20  and  0.80.  If  there  is  a  waiting  line  at  the  LIP  point 
or  if  the  IP-LIP  air  leg  is  at  full  capacity  the  touch  and 
go  aircraft  will  either  continue  to  LIP  or  re-enter  the  air 
traffic  pattern  at  the  EP1  entry  point  with  associated 
probabilities  of  0.70  and  0.30. 

It  is  also  possible  that  the  flow  of  aircraft  in  the  air 
traffic  landing  pattern  can  be  interrupted  by  emergency 
events.  Emergency  events  such  as  engine  failure,  oil 
pressure  failure,  low  fuel,  or  landing  gear  failure,  are 
common  occurrences  that  require  the  aircraft  with  the 
emergency  to  land  as  soon  as  possible.  When  such  an 
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emergency  event  occurs,  the  aircraft  with  the  emergency  pre¬ 
empts  all  other  aircraft  in  the  landing  pattern.  During  the 
emergency  the  nonemergency  aircraft  in  the  landing  pattern 
upon  reaching  the  IP,  LIP,  or  Break  air  leg  points  leave  the 
pattern  and  return  to  the  entrance  point  EP1  and  join  the 
orbiting  waiting  line.  If  a  nonemergency  aircraft  is  on  the 
Base-Final  air  leg  when  the  emergency  situation  arises  the 
aircraft  continues  and  makes  a  full  stop  landing. 

During  the  emergency  the  aircraft  present  in  the 
mission  areas  do  not  leave  their  areas  but  after  completing 

their  scheduled  activities  orbit  at  a  lower  altitude  until 
the  emergency  ends  to  avoid  interfering  in  the  local  traffic 
pattern.  Aircraft  returning  to  the  air  traffic  pattern  will 
not  enter  the  pattern  but  will  orbit  at  the  entrance  point 
until  the  emergency  ends. 

If  an  aircraft  take  off  is  scheduled  during  an  emergency 
event  the  take-off  is  delayed  until  the  emergency  event  is 
completed. 

C.  INPUT  AND  OUTPUT 

All  input  data  for  the  model  as  well  as  suggestions 
for  model  structure  was  provided  by  the  instructor 
pilots  and  from  the  control  tower  personnel  at  the 
Kalamata  Air  Force  Base.  The  general  inputs  to  the  model 
consisted  of  individual  aircraft  characteristics  and 
performance  data,  alternate  take  off  schedules,  mission 
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area  constraints,  time  distributions  for  assignments  in  the 
mission  areas,  and  the  time  distances  between  the 
reference  points  for  the  air  legs  in  the  air  traffic 
pattern . 

The  model  outputs  contain  information  pertaining  to  the 
aircraft  of  each  squadron,  the  air  traffic  pattern  and 
mission  area  utilization.  The  following  output  is 
available. 

1.  Aircraft  total  flight  time  distribution. 

2.  Number  of  aircraft  by  squadron  in  the  east  and 
the  west  mission  areas. 

3.  Number  of  aircraft  take-offs. 

4.  Number  of  entries  in  each  mission  area. 

5.  Number  of  entries  in  each  leg  of  the  air 
traffic  pattern. 

6.  Maximum  number  of  aircraft  waiting  in  each 
queue. 

7.  Average  wait  time  for  each  queue. 

8.  Waiting  time  distributions  for  the  following 

queues : 

a.  Air  traffic  controller. 

b.  Entry  points,  EP1  and  EP2. 

c.  Initial  point,  IP. 

d.  Low  Initial  point,  LIP. 
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III.  GPSS  APPLICATION 


A.  BRIEF  GPSS  DESCRIPTION 

GPSS  is  the  General  Purpose  Simulation  System  language 
developed  by  IBM  for  modeling  and  simulating  queuing 
systems.  GPSS  was  used  to  model  the  Kalamata  Air  Force  Base 
air  traffic  operations.  The  GPSS  program  is  included  in 
Appendix  A. 

GPSS  uses  the  process  interaction  approach  for  modeling 
in  which  the  model  entities  are  either  temporary  or 
permanent.  The  temporary  entities  are  called 
transactions  and  the  permanent  entities  are  called 
facilities  and  storages.  The  transactions  represent  the 
calling  population  and  the  facilities  and  storages 
represent  the  service  centers.  Transactions  interact  with 
other  transactions  and  with  the  facilities  and 
storages.  In  the  Kalamata  Air  Force  Base  model  the 
calling  populations  of  aircraft  are  represented  by 
transactions  and  the  mission  areas  and  the  air  traffic 
landing  pattern  segments  are  represented  by  facilities 
and  storages. 

The  modeling  and  programming  approach  in  GPSS  is  to 
define  a  set  of  programming  statements  called  blocks 
that  represent  the  entrance  and  flow  of  the 
transactions  into  the  queuing  system  composed  of  the 
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can  be  many 


facilities  and  storages.  There 

transactions  simultaneously  moving  through  the  blocks. 
At  any  point  in  time  each  transaction  is  positioned  at  a 
block  and  most  blocks  can  hold  many  transactions 
simultaneously.  The  transfer  of  a  transaction  from  one 
block  to  another  occurs  instantaneously  at  a  specific  time 
or  when  some  change  of  system  condition  occurs.  Time  in 
the  GPSS  model  is  managed  by  the  next  event  sequence 
with  the  simulation  clock  changing  at  nonuniform  discrete 
time  points  when  the  state  of  the  system  changes. 
Transactions  continue  to  move  through  the  system  until  they 
either  encounter  a  waiting  line  or  service  time  delay. 

In  GPSS  simulated  clock  time  is  an  integer  value 
whose  scale  value  is  chosen  by  the  programmer.  The  unit 
of  time  is  not  specifically  stated  but  is  implied  by 
providing  all  times  in  the  same  units.  In  the  Kalamata  Air 
Force  Base  model  the  unit  of  time  used  is  the  second. 

B.  MODEL  STRUCTURE  IN  GPSS 

A  brief  description  of  some  of  the  important 
programming  blocks  and  storage  areas  used  in  the  GPSS 
Kalamata  Air  Force  Base  model  are  contained  in  this 
section. 

GENERATE  and  TERMINATE  blocks:  Transactions  are  created 
and  enter  the  system  at  one  or  more  GENERATE  blocks  and 
are  removed  from  the  simulation  at  TERMINATE  blocks. 


The  time  and  frequency  with  which  transactions  enter  the 
system  are  controlled  by  the  GENERATE  block.  In  the 
Kalamata  Air  Force  Base  model  GENERATE  blocks  are  used  for 
the  entry  of  training  aircraft  from  the  squadrons  for 
take  off  assignments  and  for  the  entry  of  occasional 
aircraft  from  other  bases  into  the  Kalamata  air  traffic 
landing  pattern. 

ADVANCE  block:  The  ADVANCE  block  will  hold 

transactions  for  a  specified  or  computed  number  of  time 
units.  The  purpose  of  the  ADVANCE  block  is  to  hold 
the  transactions  in  service.  In  the  Kalamata  Air  Force 
Base  model  ADVANCE  blocks  are  used  to  simulate  the  time 
delays  associated  with  take  off  delays,  training  mission 
areas  and  transit  from  point  to  point  in  the  air  leg 
segments  of  the  air  traffic  landing  pattern. 

TEST  block:  The  TEST  block  is  used  to  manage  or  transfer 
transactions  based  upon  the  test  conditions.  In  the 
Kalamata  Air  Force  Base  model  the  TEST  block  is  used  to 
prevent  aircraft  from  entering  the  system  after  the  daily 
training  period  and  to  assure  that  the  aircraft  in  the 
system  are  correctly  processed  in  order  to  complete  all 
landings  after  the  daily  training  period  ends. 

GATE  block:  The  GATE  block  is  used  as  a  gate  to 

interrupt  the  flow  of  transactions  depending  upon 
conditions  that  set  the  gate  to  "open"  or  "closed" .  In  the 
Kalamata  Air  Force  Base  model  GATE  blocks  are  used  to 
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prevent  aircraft  from  continuing  in  the  air  traffic  landing 
pattern  or  the  take-off  queue  during  an  emergency  event. 

SELECT  block:  The  SELECT  block  is  used  to  direct  the 
flow  of  transactions.  In  the  Kalamata  Air  Force  Base 
model  SELECT  blocks  are  used  to  assign  the  aircraft  to  the 
training  mission  areas  after  take-off. 

JOIN,  REMOVE,  COUNT  and  MARK  blocks:  The  JOIN,  REMOVE, 
COUNT  and  MARK  blocks  are  used  to  collect,  remove, 
count  and  identify  transactions  in  the  queue.  In  the 
Kalamata  Air  Force  Base  model  if  the  training  mission  areas 
are  occupied  the  JOIN,  REMOVE,  COUNT  and  MARK  blocks 
are  used  to  determine  current  aircraft  assignments  based 
upon  the  shortest  waiting  lines  for  the  areas. 

FACILITIES  and  STORAGE  areas:  GPSS  FACILITIES  and 
STORAGE  areas  are  used  to  collect  and  hold  transactions 
for  time  delays  that  can  be  associated  with  service  or 
performance  of  the  transactions.  A  GPSS  FACILITY  can 
hold  one  transaction.  A  GPSS  STORAGE  area  can  hold  more 
then  one  transaction.  In  the  Kalamata  Air  Force  Base 
model  eight  FACILITIES  model  the  eight  mission  areas,  one 
FACILITY  models  emergency  aircraft,  one  FACILITY  models 
the  air  traffic  controller,  four  facilities  model 
aircraft  synchronization,  and  six  STORAGE  AREAS  have  been 
used  to  model  the  air  leg  segments  of  the  air  traffic 
landing  pattern  as  described  in  the  following  page. 
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C.  GATHERING  STATISTICS  WITH  GPSS 

GPSS  automatically  records  data  and  collects  queue 
statistics  for  transactions  that  pass  through  a  storage 
area.  In  addition  to  the  previously  described  GPSS 
programming  blocks  that  manage  the  flow  of  transactions 
there  are  several  block  types  that  are  specifically  designed 
to  gather  statistics  on  transactions.  These  blocks  and 
their  application  are  described  in  this  section. 

QUEUE  and  DEPART  statistic  blocks:  QUEUE  and  DEPART 
blocks  are  used  to  identify  specific  data  collection 
points  in  the  queue. 

Data  on  transactions  that  move  through  the  queue  and 
enter  and  leave  associated  QUEUE  and  DEPART  blocks  is 
collected  as  separate  sets  of  queue  statistics.  In  the 
Kalamata  Air  Force  Base  model  queue  statistics  were 
accumulated  over  the  entire  period  of  the  simulation,  in 
each  queue  of  the  air  traffic  pattern  and  in  each  queue  of 
the  mission  areas.  The  accumulated  statistics  for  each  of 
the  model  queues  are  identified  in  the  following  list. 

Queue  Statistics 

Maximum  number  of  aircraft. 

Average  number  of  aircraft 

Standard  deviation  of  the  number  of  aircraft. 

Average  waiting  time. 


Number  of  entries. 


Number  of  entries  that  did  not  wait  for  entry. 

Average  wait  time  to  enter  the  queue. 

TABULATE,  TABLE  and  QTABLE  statistic  blocks:  The 
TABULATE,  TABLE  and  QTABLE  blocks  are  used  to  collect 
data  for  frequency  and  cumulative  frequency  tables. 
The  TABLE  and  QTABLE  blocks  define  the  transaction 
characteristics  that  are  to  be  counted  and  the  range  of 
the  frequency  tables.  Data  on  characteristics  of  the 
transactions  that  pass  through  a  TABULATE  block  are 
automatically  collected.  The  TABULATE  block  may  be  used 
anywhere  in  the  GPSS  program  and  can  collect  data  on  all 
of  the  transactions  computed  characteristics.  The 
TABULATE,  TABLE  and  QTABLE  blocks  are  used  in  the  Kalamata 
Air  Force  Base  model  to  tabulate  the  total  flight  time 
distribution  for  the  two  types  of  aircraft  and  waiting  time 
distributions  for  all  of  the  model  queues. 


IV.  MODEL  VALIDATION 


The  Kalamata  Air  Force  Base  model  was  validated  by 
comparing  model  output  to  historical  data  for  a  specific 
take  off  schedule  for  the  squadrons  of  aircraft.  Historical 
data  was  made  available  for  each  aircraft  for  a  three  day 
period.  The  historical  data  is  included  in  Appendix  C. 

For  each  aircraft  the  historical  data  consisted  of  the 
inter  arrival  times  for  take  off,  the  time  spent  waiting  in 
the  air  traffic  controller  queue,  the  time  spent  waiting  to 
enter  area  3  and  the  total  flight  time.  The  historical  data 
was  accumulated  and  averaged  for  the  three  days  of  base 
activity.  Using  the  historical  interarrival  rates  as  model 
inputs  the  model  generated  waiting  times  and  total  flight 
times  were  compared  to  the  historical  data. 

The  historical  interarrival  rate  appeared  to  be 
nonhomogeneous  on  a  daily  basis.  Therefore  the  operations 
day  was  broken  into  time  periods  for  which  the  interarrival 
rates  were  homogeneous.  These  time  periods  were: 


1) 

07:30 

-  08:40 

5) 

11:40  - 

14:00 

2) 

08:40 

-  09:10 

6) 

14:00  - 

14:30 

3) 

09:10 

-  11:10 

7) 

14:30  - 

16:20 

4) 

11:10 

-  11:40 

8) 

16:20  - 

17:30 
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For  each  time  period  the  interarrival  times  appeared  to 
have  an  exponential  distribution.  This  hypothesis  was 
tested  using  the  Kolmogorov-Smirnov  test  and  was  not 
rejected.  The  hypothesis  test  results  are  included  in 
appendix  D.  With  this  base  information  the  model  was  run 
for  a  period  of  eight  days  using  the  hypothesized 
exponential  distributions  for  the  take  off  interarrival 
times.  The  validation  comparison  tests  follow. 

Figure  4.1  displays  a  plot  of  the  comparison  of  the 
actual  wait  time  with  model  generated  wait  time  for  entry  to 
area  3,  and  Figure  4.1a  displays  the  regression  line  for  the 
regression  model: 

Simulated  Data*A+B* ( Actual  Data). 

The  slope  of  the  regression  line  is  .97  with  a  standard 
deviation  of  .027.  The  hypothesis  B«1  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .16. 

Figure  4.2  displays  a  plot  for  the  comparison  of  actual 
wait  time  with  model  generated  wait  time  in  the  air  traffic 
controller  queue  and  Figure  4.2a  displays  the  regression 
line  for  the  regression  model: 

Simulated  Data-A+B*  (Actual  ~>ata) . 

The  slope  of  the  regression  line  is  .98  with  a  standard 
deviation  of  .043.  The  hypothesis  B*1  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .43. 
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Figure  4.3  displays  a  plot  of  the  comparison  of  actual 
total  flight  time  with  model  generated  total  flight  time 
and  Figure  4.3a  displays  the  regression  line  for  the 
regression  model: 

Simulated  Data=A+B*( Actual  Data). 

The  slope  of  the  regression  line  is  .96  with  a  standard 
deviation  of  .019.  The  hypothesis  B-l  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .09. 

These  results  provide  a  validation  of  the  model  and  show 
that  the  model  data  compares  favorably  with  the  historical 
data  for  the  specified  take  off  schedules. 


degression  of  simulated  on  actual  data 
waiting  tines  in  area  3 


Figure  4.1a  Regression  of  Sinulated  Wait  Tine 
on  Actual  Walt  Time  in  Area  3 
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Figure  4.2  Comparison  of  Actual  Wait  Time 
with  Simulated  Wait  Time  at  Line  Up 


ACTUAL  vs  SIMULATED  TOTAL  FLIGHT  Tiff 
Aircraft  T37 


(XIOO) 

IAITMC  TIME 


Figure  4.3  Comparison  of  Actual  Flight  Time 
with  Simulated  Flight  Time 
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Figure  4.3a  Regression  of  Simulated  Flight  Time 
on  Actual  Flight  Time 
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V. 


This  section  contains  a  description  of  four  basic  and 
two  additional  take-off  schedules  for  the  squadron  of 
training  aircraft  for  Kalamata  Air  Force  Base  that  were  run 
in  the  model  and  a  comparison  analysis  of  the  resulting  air 
traffic  controller  congestion  in  order  to  determine  a 
reasonable  efficient  take-off  schedule.  In  order  to  run 
each  take-off  schedule  in  the  model,  the  GPSS  program  had  to 
be  modified.  These  program  modifications  are  included  in 
Appendix  B. 

A.  SCHEDULES 

Schedule  1:  This  schedule  consists  of  a  forty  minute 
take-off  cycle.  The  cycle  consists  of  two  consecutive 
fifteen  minute  take-off  periods  followed  by  ten  minutes  of 
no  take-off  activity.  The  T-37  aircraft  squadron  assigns 
five  aircraft  for  take-off  in  the  first  fifteen  minute 
period.  The  T-2  aircraft  squadron  assigns  five  aircraft  for 
take-off  in  the  second  fifteen  minute  period.  This  cycle  is 
repeated  until  all  the  training  aircraft  are  scheduled  for 
take-off.  Aircraft  that  can  not  take-off  in  their  assigned 
period  are  recycled  for  later  take-off.  Aircraft  taking-off 
have  priority  use  of  the  runway  over  normally  landing 
aircraft.  This  schedule  was  proposed  for  analysis  by  the 
two  squadrons  at  Kalamata  Air  Force  Base. 
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Schedule  1A:  This  schedule  is  derived  from  and  is 
identical  to  Schedule  1  except  that  the  number  of  aircraft 
scheduled  for  take-off  in  each  fifteen  minute  period  is  four 
instead  of  five.  This  schedule  was  suggested  as  a  result  of 
reviewing  the  model  output  data  for  Schedule  1.  It  was 
thought  that  this  change  would  decrease  the  number  of 
training  flights  while  also  decreasing  the  air  traffic 
congestion  and  waiting  times. 

Schedule  IB:  This  schedule  is  derived  from  and  is 
identical  to  Schedule  1  with  the  take-off  cycle  period 
extended  to  sixty  minutes.  This  schedule  was  also  suggested 
as  a  result  of  reviewing  the  model  output  data  for  Schedule 
1.  It  was  thought  that  this  change  would  also  decrease  the 
number  of  training  flights  while  decreasing  the  air  traffic 
congestion  and  waiting  times. 

Schedule  2:  This  schedule  consists  of  an  eight-minute 
take-off  cycle.  The  cycle  consists  of  two  consecutive 
four-minute  take-off  periods.  The  T-37  aircraft  squadron 
assigns  one  aircraft  for  take-off  in  the  first  four-minute 
period.  The  T-2  aircraft  squadron  assigns  one  aircraft  for 
take-off  in  the  second  four-minute  period.  This  cycle  is 
repeated  until  all  the  training  aircraft  are  scheduled  for 
take-off.  Aircraft  that  can  not  take-off  in  their  assigned 
four-minute  period  are  recycled  for  later  take-off.  This 
schedule  was  also  suggested  as  a  result  of  reviewing  the 
model  output  data  for  Schedule  1.  It  was  thought  that  this 


model  output  data  for  Schedule  1.  It  was  thought  that  this 
change  would  increase  the  number  of  training  flights  while 
decreasing  the  air  traffic  congestion  and  waiting  times. 

Schedule  2A:  This  schedule  is  derived  from  Schedule  2  and 
is  identical  to  Schedule  2  except  that  the  take-off  cycle  is 
extended  to  ten  minutes  with  two  five-minute  take-  off 
periods.  This  schedule  was  suggested  as  a  result  of 
reviewing  the  model  output  data  for  Schedule  2. 

Schedule  3:  This  schedule  is  a  reproduction  of  the 

schedule  currently  in  use  at  Kalamata  Air  Force  Base.  The 
schedule  contains  no  structure  and  take-off  times  are 
scheduled  at  random.  The  distributions  of  current 
interarrival  times  for  take-off  times  were  analyzed  and 
found  to  be  exponential  for  different  periods  during  the 
day.  This  analysis  is  contained  in  Appendix  C  and  was  also 
used  for  model  validation. 

B.  RESULTS 

The  model  output  data  for  each  schedule  are  contained  in 
Appendix  E.  The  measures  of  effectiveness  used  in  comparing 
schedules  were  performance  and  efficiency/safety. 
Performance  is  measured  by  the  average  number  of  training 
aircraft  scheduled.  Efficiency  is  measured  by  the  average 
waiting  times  in  the  mission  area  queues.  Safety  is 
directly  related  to  efficiency  in  that  the  smaller  waiting 
times  mean  less  hazardous  flying  conditions.  The  empirical 
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distributions  for  both  of  these  measures  for  each  schedule 
are  contained  in  Appendix  E. 

The  summary  results  for  Schedule  1  are  of  interest  because 
this  schedule  was  suggested  by  the  training  squadrons  at 
Kalamata  Air  Force  Base.  This  schedule,  by  comparison  with 
the  other  schedule  results,  does  not  provide  the  highest 
values  for  performance  and  efficiency/safety.  It  was 
because  of  these  results  that  the  other  schedules  were 
derived  from  Schedule  1  by  sensitivity  analysis. 

The  summary  results  for  Schedule  3  are  of  interest 
because  this  schedule  is  a  reproduction  of  the  schedule 
currently  in  use  at  Kalamata  Air  Force  Base.  This  schedule 
by  comparison  with  the  other  schedule  results,  also  does  not 
provide  the  best  values  for  performance  and 
efficiency/safety. 

It  was  found  that  by  applying  Schedule  1A  and 
comparing  the  results  to  Schedule  1  that  performance 
decreased  but  that  efficiency/safety  improved  drastically. 
Figure  5.1  compares  the  cumulative  mission  area  waiting  time 
distributions  and  Figure  5.2  compares  the  cumulative  LIP 
point  waiting  time  distributions  for  Schedules  1  and  1A. 

Schedule  2  results  were  more  efficient  than  Schedule  1 
results.  However,  by  applying  Schedule  2A  mission  area 
waiting  times  decreased  dramatically  and  provided  the  most 
efficiency/safety  measure  for  all  schedules  examined. 
Figure  5.3  compares  the  cumulative  mission  area  waiting  time 
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distributions  and  figure  5.4  compares  the  cumulative  LIP 
point  waiting  time  distributions  for  schedules  2  and  2A. 
Figure  5.5  compares  the  cumulative  mission  area  waiting  time 
distributions  for  schedules  1A,  2A  and  3.  It  is  obvious 
that  Schedule  2A  is  preferable  to  Schedule  1. 

Schedule  3  results  were  more  efficient  than  Schedule  1 
results.  The  performance  measure  for  Schedule  3  however  is 
the  lowest  of  all  the  schedules  examined.  Figure  5.6 
compares  the  cumulative  LIP  point  waiting  time  distributions 
for  schedules  1A,  IB,  2A  and  3. 

C.  ANALYSIS 

From  the  above  summary  results  the  preferred  schedules 
appear  to  be  Schedules  1A  and  2A.  An  analysis  of  variance 
was  performed,  using  the  function  "ANOVA"  from  the  OA3660 
APL  WORKSPACE,  to  test  the  hypothesis  that  the  mean  mission 
area  waiting  time  differences  for  schedules  1,  1A,  2,  2A  are 
not  significant.  The  analysis  of  variance  results  are 
contained  in  Table  5.1.  These  results  show  that  the  null 
hypothesis  of  no  significant  differences  between  mean 
mission  area  waiting  times  is  rejected  at  significance 
levels  greater  than  .995. 

The  Sum  of  Squares  from  the  previous  analysis  was  broken 
into  three  components  in  order  to  test  for  individual 
effects  rather  than  just  a  schedule  effect  using  individual 
degrees  of  freedom.  The  results  of  this  analysis  are 
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contained  in  the  tables  5.2  and  5.2a.  These  results  show 
that  the  hypothesis  of  no  significant  difference  between  the 
compared  mean  mission  area  waiting  times  for  the  selected 
schedules  is  rejected  for  each  comparison  at  significance 
levels  greater  than  .975. 


TABLE  5.1 

ANALYSIS  OF  VARIANCE  RESULTS 


SOURCE  DF 

SCHEDULE  3 

ERROR  4 

TOTAL  7 

R-SOUARE  =  0.954 


OVERALL  MEAN  = 
TREATMENT  EFFECTS 
-68.5  24.02 

68.5  -24.02 


ANOVA  TABLE 

SS 

232864.74 

1X181.05 

244045.78 

200.68 

-87.56266.68* 
0.435  -17.92 
-0.435  17.92 


MS  F 

77621.59  27.77 

2795.27 


*.50  14.39 
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TABLE  5.2 


ANOVA  WITH  INDIVIDUAL  DEGREES  OF  FREEDOM 
SCHEDULE  1  SCHEDULE  1A  SCHEDULE  2A  SCHEDULE  2 


44.62 

181.62 

49.138 

443.34 

7.62 

6.75 

197.15 

232.99 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

-1 

1 

1 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

D  =  Sbij2  Z  =  S\iijXi  wi2=Zi2/Di 


8 

— 

— 

8 

716.45 

64162.575 

4 

-708.48 

125485.970 

4 

415.77 

43216.173 

2 

-137.00 

9384.500 

2 

48.04 

1153.920 

2 

0.87 

0.378 

2 

1284.50 

642.250 

TABLE  5.2a 


ANALYSIS  OF  VARIANCE  RESULTS 
USING  INDIVIDUAL  DEGREES  OF  FREEDOM 


SCHEDULES 

SS 

DF 

MS 

F 

( 1 , 1A )vs( 2 , 2A ) 

64162.575 

1 

64162.575 

22.954 

1  vs  1A 

125485.970 

1 

125485.970 

44.892 

2  vs  2A 

43216.173 

1 

43216.173 

15.460 

RESIDUAL 

11181.048 

4 

2795.262 

_ 

TOTAL 

244045.760 

| 

- - - 
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Figure  5.1  Cumulative  Waiting  Tim'?  Distribution 
in  the  Mission  Areas 


Schedules  1  and  1A 


SCHEDULES  2  AND  2A  vs  WHINS  TIME 


0  3  6  9  12  15 

(X100) 

WIT!  NS  M 


Figure  5.3  Cumulative  Waiting  Time  Distribution 
in  the  Mission  Areas 
Schedules  2  and  2A 
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Figure  5.4  Cumulative  LIP  Waiting  Times 
Schedule  2  and  2A 


39 


Figure  5.5  Cumulative  Waiting  Time  Distribution 
in  the  Mission  Areas 
Schedules  1A  2A  and  3 
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Figure  5.6  Cumulative  LIP  Waiting  Times 
Schedules  1A  IB  2A  and  3 
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VI .  CONCLUSIONS  -  SUMMARY 

A.  CONCLUSIONS 

If  an  added  emphasis  is  to  be  placed  upon  the 
efficiency/safety  factor  in  scheduling  the  aircraft  then 
Schedule  2A  is  the  preferred  schedule.  If  Schedule  2A 
conflicts  with  local  base  operations  due  to  other  operating 
constraints  then  Schedule  1A  is  the  next  best  schedule.  The 
difference  in  mean  mission  area  waiting  times  for  these  two 
schedules  is  made  more  significant  when  it  is  also 
realized  that  landing  aircraft  have  runway  priority  in 
schedule  2A  and  do  not  in  Schedule  1A.  Figure  6.1  compares 
the  cumulative  mission  area  waiting  time  distributions  and 
figure  6.2  compares  the  cumulative  LIP  point  waiting  time 
distributions  for  schedules  1A  and  2A. 

If  added  emphasis  is  to  be  placed  upon  the  performance 
factor  in  scheduling  the  aircraft  then  Schedule  2  is  the 
preferred  schedule.  However,  the  increase  of  twenty  four 
scheduled  aircraft  is  at  the  expense  of  more  than  a  ten  fold 
increase  in  mean  mission  area  waiting  times  that  contributes 
to  air  traffic  congestion  and  pilot  fatigue  and  should  be 
avoided  if  possible.  Again,  if  Schedule  2  conflicts  with 
local  base  operations  due  to  other  operating  constraints 
then  Schedule  1  is  the  next  best  schedule.  It  is 
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recommended  that  the  efficiency/safety  factor  be  the 
deciding  factor  in  selecting  schedules. 

Table  6.1  contains  a  summary  comparison  between  the 
different  schedules.  In  this  table  the  first  column 
contains  the  average  number  of  flights,  the  second  column 
contains  the  percentage  of  the  aircraft  that  waited  in  the 
mission  area  queues,  the  third  column  contains  the 
percentage  of  the  aircraft  that  waited  more  than  180  seconds 
in  the  mission  area  queues,  the  fourth  column  contains  the 
percentage  of  the  aircraft  that  waited  more  than  30  seconds 
at  LIP  and  the  fifth  column  contains  the  average  conditional 
waiting  time,  in  seconds,  in  the  mission  area  queues. 


TABLE  6.1 

SUMMARY  COMPARISON  BETWEEN  SCHEDULES 


#  of  A/C 
scheduled 

%  of  A/c 
waiting 
in 

m.  areas 

%  of  A/C 
waiting 
>18Q"in 
m. areas 

%  of  A/C 
waiting 
>30"in 
LIP 

Avg . Conditional 
waiting  time 
in  seconds 

Sch. 2 A 

106 

2 

0 

0.8 

11 

Sch.lA 

106 

7 

5.6 

1.6 

45 

Sch.  2 

130 

26 

18 

2.5 

203 

Sch. IB 

89 

7 

7 

2 

230 

Sch.  3 

99 

14 

13 

3 

368 

Sch.  1 

130 

26 

17 

3.5 

496 
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Figure  6.3  compares  the  cumulative  mission  area  waiting 
time  distributions  and  figure  6.4  compares  the  cumulative 
LIP  point  waiting  time  distributions  for  schedules  1  and  2. 


Figure  6.1  Cumulative  Waiting  Time  Distribution 
in  the  Mission  Areas 


Schedules  1A  and  2A 


Figure  6 . 


3  Cumulative  Waiting  Time  Distribution 
in  the  Mission  Areas 
Schedules  1  and  2 
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B.  SUMMARY 


This  thesis  presents  an  IBM-PC  GPSS  model  for  studying 
aircraft  take-off  schedules  at  Kalamata  Air  Force  Base  in 
Greece.  The  model  is  specific  to  Kalamata  Air  Force  Base 
because  of  the  structure  and  sequence  of  the  queues  used  in 
modeling  the  air  traffic  control  points  and  the  training 
mission  areas.  The  GPSS  program  can  be  transferred  to  any 
IBM-PC  compatible  computer  that  will  run  GPSS.  The  model 
can  therefore  be  used  at  the  Kalamata  Air  Force  Base  to 
continue  to  examine  other  aircraft  take-off  schedules  end 
to  help  reducing  fuel  consumption. 

Six  aircraft  take-off  schedules  were  examined  and  a 
comparison  of  results  was  based  upon  factors  of  performance 
and  efficiency/safety.  The  take-off  schedule  currently  used 
at  Kalamata  Air  Force  Base  is  shown  to  be  a  poor  performer 
with  a  low  efficiency/safety  factor.  The  take-off  schedule 
proposed  for  use  by  the  training  squadron  personnel  at 
Kalamata  Air  Force  Base  is  shown  to  have  the  highest 
performance  factor  and  the  worst  efficiency/safety  factor. 

The  specific  schedule  of  ten  minutes  of  a  ten  minute  take¬ 
off  cycle  consisting  of  two  consecutive  five  minute  take-off 
periods  from  each  squadron  with  runway  priority  given  to 
landing  aircraft  is  shown  to  be  the  best  schedule  based  upon 
the  recommended  emphasis  placed  upon  efficiency/safety 
rather  than  on  numerical  performance. 
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APPENDIX  A 


This  appendix  contains  the  GPSS  program  for  the  Kaiamata 
Air  Force  Base  model . 


10  ft********************************************************************* 


20  *  * 

30  *  AIR-TRAFFIC  CDNTRQL  SIMULATION  * 

40  *  * 

50  *  * 


60  ****#****#*********+#*♦***+*******#******«-*#****###**#*#************#* 

70  * 

eo  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXP  FUNCTION  RN1,C24 

0,0/. 1,. 104/. 2,. 222/. 3,. 355/. 4.. 509/. 5.. 69/. 6.. 915/. 7, 1.2/. 75, 1.38 
.8, 1.6/. 34, 1.33/. S3, 2. lZ/. 9.2. 3/. 92, 2. 52/. 94, 2. 81/. 95,2. 99/. 9o, 3. 2 
.  97, 3.  5/. 98, 3. 9/. 99, 4.  6/.  995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120  * 

130  * 

140  MEAN  FUNCTION  Cl.DS 

3600 , 360 . 1 1 /540C ,221.60/1 3200 ,488.62 

1500l, 242. 63/23400, 557. 83/25200, 124.58 

3 1 300 , 494.60/36600 ,221.30 

150  * 


160 

* 

1  30 

GENERATE 

I”1 

a 

; CREATE  A  SINGLE  TRANSACTION 

130 

SWITCH!  LOGIC  S 

; SET  LOGIC  SWITCh  1 

190 

ADVANCE 

300 

; T-37  TAKE  OFF 

200 

LOGIC  r 

1 

{RESET  LOGIC  SWICH  1 

r*  * 

ADVANCE 

2100 

{NONE  T-37  TAKES  OFF 

220 

TRANSFER 

,SWITCH1 

230 

GENERATE 

v 900,1 

240 

SWITCH2  LOGIC  S 

; SET  LOGIC  SWICH  2 

250 

ADVANCE 

300 

•T-2  TAKE  OFF 

260 

LOGIC  R 

; RESET  LOGIC  SWICH  2 

270 

ADVANCE 

2100 

{NONE  T-2  TAKES  OFF 

280 

TRANSFER 

.SWITCH2 

290 

GENERATE 

9000, FN^E 

XP 

{CREATE  EMERGENCY  EVENTS 

300 

SEIZE 

DANGE 

310 

ADVANCE 

600.180 

{EMERGENCY  HOLDS 

320 

RELEASE 

DANGE 

330 

340 

350 

TERMINATE 

{EMERGENCY  TERMINATES 

GENERATE 

9100.FNTEXP 

; A/C  FROM  OTHER  BASES 

360 

ASSIGN 

1  O 

*  t  *. 

{ASSIGN  TO  A  PARAMETER 

370 

MARK 

4 

{MARK  THE  TIME 

380 

TRANSFER 

,EPI 

390 

• 

400 

*  CREATE  AIRCRAFTS 

FROM  TWO 

DIFFERENT 

SQUADROMS  ON  BASE 

410 

* 

420 

GENERATE 

60.10 

{CREATE  T-37 'S 

430 

TEST  L 

N$GNkT , 1 , 
1  ,PISTA 

PI  ST  A 

;THE  SIMULATION  STOPS  AFTER 

440 

SATE  LS 

450 

ASSIGN 

1,1 

{ASSIGN  TO  A  PARAMETER 

460 

TRANSFER 

,  DOWN 

470 

* 

430 

GENERATE 

60.10 
NfGNRT, 1 , 

{CREATE  T-2 'S 

490 

TEST  L 

PI  ST  A 

500 

GATE  LS 

2,PISTA 
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5:0 

ASSIGN 

1,0 

; ASSIGN  ~D  A  PARAN.ETEF 

520  DOWN 

QUEUE 

CNTR 

; WAIT  0UT  OR  THE  RUNWAY 

C7-?  i' 

GATE  FV 

DANSE 

; DON'T  MOVE  IF  EMERGENCY  HCOE 

540  ¥■ 

GATE  RE 

CAP4 

; CHECK  FOP  A\C  ON  EASE  LEZ 

550 

SEIZE 

CNTF 

;  CAPTURE  THE  COtTRCLLE-. 

560 

DEPART 

CNTR 

; GOING  FOR  LINE  UP 

570  LINUF 

ADVANCE 

70,20 

CNTR 

; LINE-UP  CHECK 

580 

RELEASE 

; TAKE-OFF 

590 

MARK 

6 

: START  FLIGHT  TIMEFOR  T-2  A/C 

600 

MARK 

7 

; START  FLIGHT  TIME  FOR  T-37  A/C 

610 

ADVANCE 

100,10 

; AFTER  TAKE  OFF  TO  DEPART  POINT 

615 

620  * 

TRANSFER 

.25,, OTHER 

{FORMATIONS,  INSTR.FL. ,  CPM 

630 

SELECT  E 

2,101 ,108,0, F.QUEUP 
BV$F I RST , i , WEST 

:FIND  EMPTY  AREA  IF  EXISTS  ONE 

640  AREA 

650  * 

TEST  E 

{CHOOSE  EAST  OR  WEST  AREA 

660  *  E 

670  * 

AST  A  R 

E  A  S 

680 

SAVEVALUE 

13+,P1 

{RECORDS  T-37 'S 

690 

SAVEVALUE 

14+,  1 

{RECORDS  TOTAL  #  OF  A/C 

700  WAIT 

QUEUE 

P2 

{WAIT  IN  THE  QUEUE 

710 

TEST  E 

P1,1,TTW0 

{A/C  TYPE  T-37 

720 

SEIZE 

P2 

{CAPTURE  THE  AREA 

730 

DEPART 

P2 

{LEAVE  THE  QUEUE 

740 

JOIN 

FARM 

{JOIN  THE  GROUP 

750 

MARK 

60 

{COUNT  THE  TIME  IN  THE  AREA 

760 

ADVANCE 

2100,300 

{DO  SCHEDULED  ACTIVITY 

770 

RELEASE 

P2 

{FREE  THE  AREA 

780 

790  * 

REMOVE 

FARM 

{OUT  OF  THE  GROUP 

800 

GATE  FV 

DANGE 

{IF  EMERGENCY  ORBIT  UNDER 

810  TOIP 

ADVANCE 

120,30 

cap! 

{GOING  TO  EPI 

820  I PE 

QUEUE 

{ENTER  THE  QUEUE 

830 

SEIZE 

DUMY1 

{SYNCHRONIZE  THE  AIRCRAFTE 

840 

GATE  FV 

DANGE 

; IF  EMERGENCY  ORBIT  THERE 

850 

ENTER 

CAPl 

{CAPTURE  THE  STORAGE 

860 

TRANSFER 

SIM, DIR, DLY 

870  DLY 

ADVANCE 

20 

880  DIR 

DEPART 

CAPl 

{LEAVE  THE  QUEUE 

890 

RELEASE 

DUMYI 

900 

ADVANCE 

120,10 

CAPl 

{GOING  TOWARDS  EP2 

910 

920  * 

LEAVE 

{FREE  THE  STORAGE 

930  BAK 

QUEUE 

CAP2 

{ENTER  THE  QUEUE 

940 

SEIZE 

DUMY2 

{SYNCHRONIZE  THE  AIRCRAFTS 

950 

GATE  FV 

DANGE 

{IF  EMERGENCY  DRBIT  THERE 

960 

ENTER 

CAP2 

{ENTER  THE  STORAGE 

970 

TRANSFER 

SIM, DRC, DLL 

980  DLL. 

ADVANCE 

10 

990  DRC 

DEPART 

CAP2 

{DEPART  THE  QUEUE 

1000 

RELEASE 

DUMY2 

1010 

ADVANCE 

70,10 

CAP2 

{GOING  TOWARDS  IP 

1020 

1030  * 

LEAVE 

{FREE  THE  STORAGE 

1040  RWAY 

QUEUE 

CAP3 

{ENTER  THE  QUEUE  AT  IP 

1050 

SEIZE 

DUMY3 

{SYNCHRONIZE  THE  AIRCRAFTS 

1060 

SATE  FV 

DANGE, EPI 

{IF  EMERGENCY  60  TO  EPI 

1070 

ENTER 

CAP3 

{ENTER  THE  STORAGE 

1080 

TRANSFER 

SIM, DIREC, DELAY 

1090  DELAY 

ADVANCE 

10 

1100  DIREC 

DEPART 

CAP3 

{DEPART  THE  QUEUE 

1110 

RELEASE 

DUMY3 

1120 

ADVANCE 

65,5 

CAP3 

; GOINS  TOWARDS  LIP 

1130 

1140  » 

LEAVE 

{FREE  THE  STORAGE 

1150  LIP 

QUEUE 

CAP4 

{ENTER  THE  QUEUE  AT  LIP 

1160 

SEIZE 

DUMY4 

{SYNCHRONIZE  THE  AIRCRAFTS 

1170 

GATE  FV 

DANGE, EPI 

{IF  EMERGENCY  GO  TO  EPI 

1180 

ENTER 

CAP4 

;  ENTER  THE  STORAGE 
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U2&  A  transfer 

;20>.  Dr_2  ADVAN^c 

D:RS  DEPART 
1-~-  RZ_EA-~ 

ADVANCE 
LEAVE 
3ATE  FV 
GATE  £NF 


1240 

1250 

1260 


1270  * 

}i§9  enter 

}290  ADVANCE 

1300  LEAVE 

1310  GATE  P u 

1320  BASE  GATE  SNF 


1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1430 
1433 
1436 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
15  30 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1840 


ENTER 

advance 
LEAVE 

GATE  FV 
GATE  FV 
TEST  NE 
TEST  E 
TEST  LE 
MORE  ADVANCE 
TES”  NE 
SAVEVALUE 
TRANSFER 
OUT  TEST  S 
TEcO  ADVANCE 
_Tt,  _  TEST  NE 

finish  test  e 
TABULATE 
„  TRANSFER 

TABULT2  TABULATE 
PISTA  TERMINATE 

CNTU  ADVANCE 
GATE  FV 
TRANSFER 
MLRt  GATE  SNF 

HERE1  TEST  E 

ADVANCE 
TRANSFER 

* 

EPI  ADVANCE 

TRANSFER 

# 

GOAR  ADVANCE 

BATE  SNE 
TRANSFER 
EPI1  TRANSFER 

* 

TTWO  SEIZE 

DEPART 
JOIN 
MARK 
ADVANCE 
RELEASE 
REMOVE 
TRANSFER 

* 

*  W  E 


1850 


WEST  SAVEVALUE 
SAVEVALUE 
QUEUE 

TEST  E 


SIh9OI  Rc.  9  DE»_  A 

CA&4 

Djrv- 

; LEAVE  THE  QUEUE 

caW 

DANCE. E- I 

CAPS, CNTU 

! GOINS  -re  BREA-'  PZIK~ 

;FREE  THE  STORAGE 

;  IF  EMERGENCY  :  GC  td  EVrs 

ji-HEwk  A/C  ON  DOWNLINE  rON*~_C_-  r 

CAPS 

45,15 

CAP5 

DANSE, FAST 

CAP6 , GOAR 

jSSLiril  SgCHR3N,2E0  ** 

5  FREE  THE  STORAGE 
» IF  EMERGENCY  :  GO  TO  ENTR  phtm- 
1  CHECK  #A/C  ON  BASELEG  DON'T  QEJEU 

CAP6 

45.15 

CAP6 

jfHTER  THE  CAPACITY 
jBASc  LEG 

{FREE  THE  CAPACITY 

DANGE.TELO 

CNTR , GOAR 

PI, 2, OUT 
Pl.l.TTT 
MP7,3900,TEL0 
10,5 

PI, 2, CNTU 

!rc  fHIRGENCY  DO  FULL  STOP 

5  CHECK  FOR  A/C  TYPE  T-37 
?  CHECK  FLIGHT  TIME  OF  T-37 
;  TOUCH  AND  GO 

20+ ,  i 

,  CNTU 

$  COUNT  THE  NUMBER  OF  LANDS 

MP4, 900, MORE 

20,5 

PI, 2, PISTA 

PI i 1 , TABULT2 

T0TAL37 

.PISTA 

T0TALT2 

5 CHECK  FLIGHT  TIME  FOR  FOF -  £  ' " 

;FULL  STOP  LANDING  ^  ~ 

! CHECK  FOR  OTHER  TYPE  A/r 
;IF  A/C  TYPE  IS  T-37  CDN’-’NIIP 
{DISTRIBUTION  OF  THE  TOTAL I ME 

{TABULATE  T-2  FLIGHT  TIME 

60,5 

DANGE.CRIT 
.20. .CLOSE 
CAP3.EP 11 
Q*CAP4,0,EPI 
120,20 
.LIP 

{CONTINUE  TO  LIP 

5 IE EMERGENCY  :  GO  TO  EN~R. PQ’\- 
» 204  CONTINUE  FOR  CLOSE  PALTER* 
.CHECK#  A/C  FROM  IP  TO  LIP  L 
{IF  A/C  WAITS  AT  LIP,  GC  TO  EPT 
{GOING  TO  LIP  u  t  - 

170,20 

.IF* 

{GOING  TC  EPI 

30,10 

CAP5.CL0SE 

.HERE 

{GOING  AROUND 

{IF  NO  A/C  IN  DOWNWING  DC  CLOSE 

. 30.HERE1 ,EPI 

{30%  GO  TO  EPI 

P2 

P2 

FARM 

60 

2400,300 

P2 

FARM 

,TOIP 

5  CAPTURE  THE  AREA 

J LEAVE  THE  QUEUE 

{JOIN  THE  GROUP 

{MARK  THE  TIME  IN  THE  AREA 

{OUT  OF  THE  GROUP 

ST  AREAS 

15+,P1 

16+,  1 

P2 

{RECORDS  T-37'S 
{RECORDS  TOTAL  #  A/C 
{ENTER  THE  QUEUE 

PI, 1,TTWB 

: CHECK  T4E  TYRE  Oc  THE  A'C 
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i860 

SEIZE 

C'2 

; A/C  T-37  CAPTURES  THE  AREA 

1870 

DE°AR" 

P2 

; LEAVE  THE  QUEUE 

133:. 

JC  IN 

FART 

; JOIN  THE  GROUP 

1 890 

MARK 

60 

•MARK  THE  TIME  IN  "HE  AREA 

<  c 

A 

ADVANCE 

2100,3c:- 

;D0  THE  SCHEDULED 

1910 

RELEASE 

P2 

;FREE  THE  AREA 

1820 

REMOVE 

FARM 

{OUT  OF  THE  GROUc' 

1930 

ADVANCE 

200,100 

; GOING  TO  EP2 

1940 

1950  * 

TRANSFER 

,  BAK 

1960  TTWB 

SEIZE 

P2 

; A/C  T-2  CAPTURES  THE  AREA 

1970 

DEPART 

P2 

;  LEAVE  THE  QUEUE 

1980 

JOIN 

FARM 

;  JOIN  THE  GROUP 

1990 

MARK 

60 

{HARK  THE  TIME  IN  THE  AREA 

2000 

ADVANCE 

2400.300 

{DO  SCHEDULED  ACTIVITY 

2010 

RELEASE 

P2 

; FREE  THE  AREA 

2020 

REMOVE 

FARM 

; OUT  OF  THE  GROUP 

2030 

ADVANCE 

200,100 

; GOING  TO  EP2 

2040 

2050  * 

TRANSFER 

,BAK 

2060  TTT 

TEST  LE 

MP^SOO ,  TELO 

{CHECK  FLIGHT  TIME  FOR  T-2  A/C 

2070 

2080  * 

TRANSFER 

;G0  FOR  FULL  STOP 

2090  FAST 

TRANSFER 

•80,EPI ,BASE 

; 207.  BD  TO  I PE,  OTHERS  CONTINUE 

2100  GRIT 

TRANSFER 

.80.FST.EPI 

CNTR.HERE 

45.15 

.BASE 

;  807.  GO  TO  IPE 

2110  CLOSE 

GATE  FV 

;  IF  NO  A/C  FOR  TAKE  OFF  CONTINUE 

2120 

ADVANCE 

{DOING  CLOSE  PATTERN 

2130 

TRANSFER 

;BASE  KEY 

2140  * 

2150  QUELP 

COUN~  GE 

'pl^'SLCt1’0 

r  lvfOj  oLLr  1 

; CHECK  THE  QUEUES  IN  ALL  AREAS 

2160 

TEST  L 

{CHECK  IF  ALL  QUEUES  ARE  OCCUPIED 

2170 

SCAN  MAX 

FARH.MP6Q, ,P2,2 

{FIND  THE  EARLIEST  OCCUPIED  AREA 

2180 

SAVEVAL'JE 

32  P  ■ 

FAfelT ALL , ,  P2 , X32 
i?l01,108,,Q 

2190 

REMOVE 

.AREA  {REMOVE  FROM  THE  GROUP 

2200 

TRANSFER 

2210  SLCT 

SELECT  MIN 

{SELECT  AREA  WITH  THE  MIN  QUEUE 

r»r>on 

A-  A,  A. 

SAVEVALUE 

32  P2 

Fa6m^ALL,,P2,X32 

2230 

REMOVE 

,AREA  {REMOVE  FROM  THE  GROUP 

2240 

TRANSFER 

2243  OTHER 

SAVEVALUE 

$o+,i 

2400,300 

{A/C  NDT  GOING  TO  MISSION  AREAS 

2246 

ADVANCE 

{PERFORM  SHEDULED  ACTIVITIES 

2248 

2249  * 

TRANSFER 

,IPE 

{BACK  TO  THE  AIRPORT 

2250  * 

2260  * 

2270  * 

M 

AKE  TABLE 

S 

2280  TOTALS 

7  TABLE 

MP7, 3600, 240, 18 

{TOTAL  TIME  DISTRIBUTION 

2290  T0TALT2  TABLE 

MP6f 3600,240, 18 

{TOTAL  TIME  DISTRIBUTION 

2300  CNTR 

QTABLE 

CNTR, 0,60, 40 

{TIME  DISTRIBUTION  IN  CNTR  QUEUE 

2310  CAP1 

OTABLE 

CAP1 ,0,30,40 

{Tift  DISTRIBUTION  IN  CAP1  QUEUE 

2320  CAP2 

QTABLE 

CAP2,0,30,40 

{TIME  DISTRIBUTION  IN  CAP2  QUEUE 

2330  CAP3 

QTABLE 

CAP3,0,30,40 

{Tilt  DISTRIBUTION  IN  CAP3  QUEUE 

2340  CAP4 

QTABLE 

CAP4,0,30,40 

{TIME  DISTRIBUTION  IN  CAP4  QUEUE 

2350  AREA1 

QTABLE 

101,0,120,32 

{TIME  DISTRIBUTION  IN  AREA  Q101 

2360  AREA2 

QTABLE 

102,0,120,32 

{TIME  DISTRIBUTION  IN  AREA  Q102 

2370  AREAS 

QTABLE 

103,0,120,32 

{TIME  DISTRIBUTION  IN  AREA  Q103 

2380  AREA4 

QTABLE 

104,0,120,32 

{TIME  DISTRIBUTION  IN  AREA  0104 
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i 


2390 

AREA! 

C7A3i_E 

10“, 0,120, 32  ; TIME  DISTkIEj  I0'« 

*  i  v 

AREA 

2400 

AREA6 

C'ABwE 

106,0,120,32  ;  TIME  DISTRIBUTION 

IN 

AREA 

f  -  ■*, 

*  V 

AREA? 

ST  ABi_E 

107,0,120,32  ; T I ME  DISTRIBUTION 

T  « 

l  1 

ARE^ 

2420 

AREAS 

STABLE 

108,0,120,32  ; TIME  DISTRIBUTION 

IK 

Ar\£_A 

247.0 

* 

2440 

FIRST 

B VARIABLE 

<  P2 ‘ E ' 1 0 1 +P2 ' E ' 1 02+P2 ‘ E ' 1 03+P2 ' E ' 1 04 ) 

2450 

* 

2460 

* 

D  E  F  I  N 

E 

THE  STORAGES 

2470 

CAP1 

STORAGE 

2 

2480 

CAP2 

STORAGE 

2 

2490 

CAP3 

STORAGE 

2 

2500 

CAP4 

STORAGE 

3 

2510 

CAPS 

STORAGE 

T 

2520 

CAP6 

STORAGE 

2 

2530 

* 

2540 

* 

TIME 

S 

E  G  M  E  N  T 

2550 

3NF:7 

GENERATE 

30600  sTIME  ARRIVES  AT  3060' 

0 

2560 

TEST  E 

Nf DOWN, N* FINISH  :WAIT  UNTIL  ALL  THE  A/C 

LAND 

2570 

TERMINATE 

1  ; SHUT  OFF  THE  RUN 
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This  appendix  contains  the  GPSS  program  modifications 
for  the  different  take  off  schedules. 

SCHEDULE  1 


10  it********************************************************************* 
20  *  ♦ 

30  *  AIR-TRAFFIC  CONTROL  SIMULATION  * 

40  *  * 

50  *  * 

60  ********************************************************************* 
70  * 

80  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXP  FUNCTION  RN1,C24 


0.0/.  1,.  104/. 2,. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6,. 915/. 7, 1.2/. 75.1. 38 
. B, 1.6/. 84. 1.83/. 83, 2. 12/. 9, 2. 3/. 92, 2. 52/. 94, 2. 81/. 95,2. 99/. 96, 3. 2 
.97, 3. 5/. 99, 3. 9/. 99,4. 6/. 995. 5. 3/. 998, 6. 2/. 999, 7/.  9998, 8 
120  * 

130  * 

140  MEAN  FUNCTION  C1.D8 

3600 ,380.11 /5400 ,221.60/1 3200 , 488. 62 

15000 , 242.63/23400 , 557 . 83/25200 ,124.58 

ADK  L  /\  t  l  f\f\  Trt 


31800, 494 . 60/36600 ,221.30 
150  * 


160  * 

170 


GENERATE 


180  SWITCH1  LOGIC  S 
190  ADVANCE 

200  LOGIC  R 

210  ADVANCE 

220  TRANSFER 

230  GENERATE 


240  SWITCH2  LOGIC  S 

250 

ADVANCE 

260 

LOGIC  R 

270 

ADVANCE 

260 

TRANSFER 

290 

GENERATE 

300 

SEIZE 

310 

ADVANCE 

320 

RELEASE 

330 

TERMINATE 

340  * 

350 

GENERATE 

360 

ASSIGN 

370 

MARK 

380 

TRANSFER 

390  * 

400  * 

CREATE  AIRCRAI 

410  * 

420 

GENERATE 

430 

TEST  L 

440 

GATE  LS 

450 

ASSIGN 

460 

TRANSFER 

470  * 

480 

GENERATE 

300 

1 

2100 
,SWITCH1 
j, 900,1 

300 

o 

2100 

•SWITCH2 

9000 .FNtEXP 

DANGE 

600,180 

DANGE 


9100,FN$EXP 

i’2 

4 

,EPI 


; CREATE  A  SINGLE  TRANSACT 
;SET  LOGIC  SWITCH  1 
;T-37  TAKE  OFF 
; RESET  LOGIC  SWICH  1 
;N0NE  T-37  TAKES  OFF 

:SET  LOGIC  SWICH  2 
; T-2  TAKE  OFF 
; RESET  LOGIC  SWICH  2 
;N0NE  T-2  TAKES  OFF 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

{A/C  FROM  OTHER  BASES 
•ASSIGN  TO  A  PARAMETER 
{MARK  THE  Tift 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


60,10 

N4GNRT. 1 ,PI5TA 
1  ,PISTA 
M 

,  DOWN 

60.10 


{CREATE  T-37'S 

{THE  SIMULATION  STOPS  AFTER 

{ASSIGN  TO  A  PARAMETER 


: CREATE  T-2'S 


SCHEDJwS  1A 


10  ******#*******************»******************************************* 


20  *  * 

30  *  AIR-TRAFFIC  CONTROL  SIMULATION  * 

50  *  * 


60  ft********************************************************************* 

70  * 

BO  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXP  FUNCTION  RN1,C24 

O.O/.l,.  104/. 2,. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6,. 915/. 7, 1. 2/. 75, 1.36 
.8, 1.6/. 84. 1.83/. 88, 2. 12/. 9,2. 3/. 92, 2. 52/. 94, 2. 81/. 95,2. 99/. 96, 3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995,5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120  * 

130  * 

140  MEAN  FUNCTION  Cl .08 

3600 ,380.11 /5400 ,221. 60/ 1 3200 ,488.62 
1 5000 , 242 . 63/23400 , 557 . 83/25200 ,124.58 
31800 , 494 . 60/36600 , 22 1 . 3C 

150  * 


160  * 

170 

GENERATE 

180  SWITCH1 

LOGIC  8 

190 

ADVANCE 

200 

LOGIC  P 

210 

ADVANCE 

220 

TRANSFER 

230 

GENERATE 

240  SWITCH2 

LOGIC  S 

250 

ADVANCE 

260 

LOGIC  R 

270 

ADVANCE 

280 

TRANSFER 

290 

GENERATE 

300 

SEIZE 

310 

ADVANCE 

320 

RELEASE 

330 

TERMINATE 

340  * 

350 

GENERATE 

360 

ASSIGN 

370 

MARK 

380 

TRANSFER 

300 

1 

210C 
,SWITCN1 
v 900,1 

300 

2 

2100 

.SWITCH2 

9000.FN4EXP 

DANGE 

600,180 

DANGE 


9100,FN*EXP 

i’2 

4 

,EPI 


: CREATE  A  SINGLE  TRANSACTION 
■SET  LOGIC  SWITCH  1 
; T-37  TAKE  D-F 
; RESET  LOGIC  SWICH  1 
•NONE  T-37  TAKES  OFF 


;SET  LOGIC  SWICH  2 
; T-2  TAKE  OFF 
•RESET  LOGIC  SWICH  2 
■NONE  T-2  TAKES  OFF 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

:  A/C  FROM  OTHER  BASES 
•ASSIGN  TO  A  PARAMETER 
{HARK  THE  TIME 


390  * 
400  * 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


410  * 

420 

430 

440 

450 


6ENERATE 
TEST  L 
GATE  LS 
ASSIGN 


75,10 

N$GNRT.1,PISTA 

1 ,PISTA 

1.1 


{CREATE  T-37 'S 

•THE  SIMULATION  STOPS  AFTER 

{ASSIGN  TO  A  PARAMETER 


460 
470  * 


480 


490 


TRANSFER  ,DOWN 

GENERATE  75.10  {CREATE  T-2'S 

TEST  L  N$SNRT,1,PISTA 


500  GATE  LS  2,PISTA 

510  ASSIGN  1,0 

520  DOWN  QUEUE  CNTR 

530  GATE  FV  DANGE 

540  *  GATE  SE  CAPfc 


9 


ASSIGN  TO  A  PARAMETER 
WAIT  OUT  OF  THE  RUNWAY 
DONT  MOVE  IF  EMERGENCY  HOLDS 
CHECK  FOR  A\C  ON  BASE  LEG 
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SCHEBU.E  IE 


1 0  ***********************************************************************  * 


20  #  * 

30  *  AIR-TRAFFIC  CONTROL  SIMULATION  * 

40  #  * 

50  *  * 


60  Ht*MtMw***m»*H*m*»***»«m*mtt*HH«tt***tM**mH»»*mM 
70  * 

80  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  # 

110  EXP  FUNCTION  RN1.C24 

0,0/. 1,. 104/. 2,. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6.. 915/. 7, 1.2/. 75, 1.38 
.8,1. 6/. 84, 1.83/. 88, 2. 12/. 9, 2. 3/. 92, 2. 52/. 94, 2. 81/. 95.2. 99/.  96, 3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 

130  * 

140  MEAN  FUNCTION  C1,D8 

3600 ,380.11/5400.221.60/1 3200 , 488 .  62 
15000,242. 63/23400,557. 83/25200 , 124. 58 
31 800 , 494 . 60/36600 ,221. 30 
150  * 

160  * 


170 

GENERATE 

,,,1 

; CREATE  A  SINGLE  TRANSACTION 

180  SWITCH1 

LOGIC  S 

1 

{SET  LOGIC  SWITCH  1 

190 

ADVANCE 

300 

;T-37  TAKE  OFF 

200 

LOGIC  R 

1 

; RESET  LOGIC  SWICH  1 

210 

ADVANCE 

3300 

; NONE  T-37  TAKES  OFF 

220 

TRANSFER 

,SWITCH1 

230 

GENERATE 

1,900,1 

240  SWITCH2 

LOGIC  S 

;SET  LOGIC  SWICH  2 

250 

ADVANCE 

300 

;T-2  TAKE  OFF 

260 

LOGIC  R 

2 

; RESET  LOGIC  SWICH  2 

270 

ADVANCE 

3300 

{NONE  T-2  TAKES  OFF 

280 

TRANSFER 

•SWITCH2 

290 

GENERATE 

9000.FN4EXP 

; CREATE  EMERGENCY  EVENTS 

300 

SEIZE 

DANGc 

310 

ADVANCE 

600.180 

; EMERGENCY  HOLDS 

320 

RELEASE 

DANGE 

330 

TERMINATE 

; EMERGENCY  TERMINATES 

340  * 

350 

GENERATE 

9100,FN$EXP 

; A/C  FROM  OTHER  BASES 

360 

ASSIGN 

1,2 

; ASSIGN  TO  A  PARAMETER 

370 

MARK 

4 

;MARK  THE  TIME 

380 

TRANSFER 

,EPI 

390  * 

400  *  CREATE  AIRCRAFTS 

FROM  TWO  DIFFERENT 

SQUADRONS  ON  BASE 

410  # 

420 

GENERATE 

60,10 

; CREATE  T-37'S 

430 

TEST  L 

N*GNRT,1,PISTA 

} THE  SIMULATION  STOPS  AFTER 

440 

GATE  LS 

1,PISTA 

450 

ASSIGN 

1,1 

{ASSIGN  TO  A  PARAMETER 

460 

TRANSFER 

,DOWN 

470  * 

480 

GENERATE 

60,10 

{CREATE  T-2 'S 

490 

TEST  L 

N$GNRT, 1 ,PISTA 

500 

SATE  LS 

2,PISTA 

510 

ASSIGN 

1,0 

; ASSIGN  TO  A  PARAMETER 

520  DOWN 

QUEUE 

CNTR 

{WAIT  OUT  OF  THE  RUNWAY 

530 

GATE  FV 

DANGE 

;DONT  MOVE  IF  EMERGENCY  HOLD 

540  * 

GATE  SE 

CAPfc 

; CHECK  FOR  A\C  ON  BASE  LEE 

SCHEDULE  2 


10  **************#***#************************#**********#******.#.******** 
20  *  * 

30  •*  AIR-TRAFFIC  CONTROL  SIMULATION  # 

40  *  * 

50  *  * 

70  * 

80  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXP  FUNCTION  RN1.C24 

0,0/.  1, . 104/. 2,. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6,. 915/. 7, 1.2/. 75, 1.38 
.8, 1.6/. 84, 1.83/. 88, 2. 12/. 9, 2. 3/. 92, 2. 52/. 94, 2. 81/. 95.2. 99/. 96, 3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998,8 
120  * 

130  * 

140  MEAN  FUNCTION  Cl,De 

3600 .380.11 /5400 ,221.60/1 3200 .488.62 

1 5000 , 242 . 63/23460 , 557 . 83 /25260 ,124.58 

31800,494.60/36600,221.30 

150  * 

160  * 


290 
300 
310 
320 
330 
340  # 

350 
360 
370 
380 
390  * 

400  *  CR£ 

410  * 

411 

412 

413 

414 

415 

416 

417  * 

420 

430 

431  * 

450 
460 
470  # 

480 

490 

491  * 

510 

520  DOWN 
530 
540 
550 

560 

561  « 

562  * 

570  LINEUP 

580 

590 

600 


GENERATE 

SEIZE 

ADVANCE 

RELEASE 

TERMINATE 

GENERATE 

ASSIGN 

MARK 

TRANSFER 


9000.FN4EXP 

DANGER 

600,180 

DANGER 


9100,FN*EXP 

{.2 

,TQEP1 


; CREATE  EMERGENCY  EVEN'S 

; EMERGENCY  HOLDS 

{EMERGENCY  TERMINATES 

;A/C  FROM  OTHER  BASES 
{ASSIGN  TO  A  PARAMETER 
{MARK  THE  TIME 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


GENERATE 

ASSIGN 

TRANSFER 

GENERATE 

ASSIGN 

TRANSFER 

GENERATE 
TEST  L 

ASSIGN 

TRANSFER 

GENERATE 
TEST  L 

ASSIGN 
QUEUE 
GATE  FV 
GATE  SE 
SEIZE 
DEPART 


ADVANCE 

RELEASE 

MARK 

MARK 


120, 10, ,4 

,5own 

120,10,60,4 

lioWN 

480,60,1500 
N4GNRT, 1 ,PISTA 

I&OWN 

480,60,1740 

N*GNRT,1,PISTA 


DANGER 

CAP6 

CNTR 

CNTR 

N  E  UP 
70,20 
CNTR 
6 
7 


{CREATE  THE  FIRST  4  T-37S 
{ASSIGN  TO  A  PARAMETER 

{CREATE  THE  FIRST  4  T-2'S 
{ASSIGN  TO  A  PARAMETER 


• TRFflTP  T-77 ' C 
{THE  SIMULATION  STOPS  AFTER 
{ALL  THE  A/C  HAVE  LANDED 
{ASSIGN  TO  A  PARAMETER 


{CREATE  T-2'S 

{THE  SIMULATION  STOPS  AFTER 
{ALL  THE  A/C  HAVE  LANDED 
{ASSIGN  TO  A  PARAMETER 
{WAIT  OUT  OF  THE  RUNWAY 
{DONT  MOVE  IF  EMERGENCY  HOLDS 
{CHECK  FOR  A\C  ON  BASE  LEG 
{CAPTURE  THE  CONTROLLER 
{GOING  FOR  LINE  UP 


{LINE-UP  CHECK 
{TAKE-OFF 

{START  FLIGHT  TIMEFOR  T-I 
{START  FLIGHT  TIME  FOF  T- 
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SCHEDULE  2A 


10  ******************************-****************************«*********** 
20  *  * 

30  *  AIR-TRAFFIC  CONTROL  SIMULATION  * 

40  *  * 

50  *  * 

60  ***********«****#***********************4HHHHHHHHHHHHHHHHt************* 

70  * 

80  * 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXP 

0,0/.  1 


FUNCTION  RN1,C24 
104/. 2,. 222/. 3,. 355/. 4,  ' 

n 


•  V/  •  i  li  1  Vt/  «4.i  t  444/  i^l  •  v^JJ/  •  T  ,.509/. 5.. 69/. 6,. 915/. 7, 1.2/. 75, 1.38 
.8, 1.6/. 84, 1.83/. 88, 2. 12/. 9, 2. 3/. 92, 2. 52/. 94, 2. 81/. 95,2. 99/. 96, 3. 2 
.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 
120  * 

130  * 

140  MEAN  FUNCTION  C1,D8 

3600 , 380 . 1 1 /5400 ,221.60/13200,488. 62 

1 5000 , 242 . 63/23400 , 557 . 83/25200 ,124.58 

31800,494.60/36600,221.30 

150  * 


160  * 

290 

300 

310 

320 

330 

340  * 

350 

360 

370 

380 

390  * 

400  * 

410  * 


GENERATE 

SEIZE 

ADVANCE 

RELEASE 

TERMINATE 

GENERATE 

ASSIGN 

MARK 

TRANSFER 


9000.FN*EXF 

DANGER 

600,180 

DANGER 


9100,FN*EXP 

V2 

, TOEP 1 


; CREATE  EMERGENCY  EVENT; 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

; A/C  FROM  OTHER  BASES 
; ASSIGN  TO  A  PARAMETER 
;MARK  THE  TIME 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


411 

GENERATE 

412 

ASSIGN 

413 

TRANSFER 

414 

GENERATE 

415 

ASSIGN 

416 

TRANSFER 

417  * 

420 

GENERATE 

430 

TEST  L 

431  # 

450 

ASSIGN 

460 

TRANSFER 

470  * 

480 

GENERATE 

490 

TEST  L 

491  * 

510 

ASSIGN 

520  DOWN 

QUEUE 

530 

GATE  FV 

540 

GATE  SE 

550 

SEIZE 

560 

DEPART 

561  * 

562  * 

570  LINEUP 

ADVANCE 

580 

RELEASE 

590 

MARK 

600 

MARK 

120, 10, ,4 

,6own 

120,10,60,4 

IftoWN 

600,60,1500 

N*GNRT,1,PISTA 


\ixm 


600,60,1800 
N4GNRT, 1 ,PISTA 

1,0 

CNTR 

DANGER 

CAP6 

CNTR 

CNTR 


L  I 


N  E 

70,20 

CNTR 

6 

7 


U  P 


; CREATE  THE  FIRST  4  T-37'S 
; ASSIGN  TO  A  PARAMETER 

; CREATE  THE  FIRST  4  T-2'S 
; ASSIGN  TO  A  PARAMETER 


•  CREATE  T-37 ' S 
| TIC  SIMULATION  STOPS  AFTER 
; ALL  THE  A/C  HAVE  LANDED 
; ASSIGN  TO  A  PARAMETER 

; CREATE  T-2'S 

;THE  SIMULATION  STOPS  AFTER 
; ALL  THE  A/C  HAVE  LANDED 
t  ASSIGN  TO  A  PARAMETER 
{WAIT  OUT  OF  THE  RUNWAY 
jDONT  MOVE  IF  EJCRGENCY  HOLDS 
{CHECK  FOR  A\C  ON  BASE  LEG 
{CAPTURE  THE  CONTROLLER 
{GOING  FOR  LINE  UP 


; LINE-UP  CHECK 
•TAKE-OFF 

{START  FLIGHT  TIMEFOR  T-2  A/C 
{START  FLIGHT  TIME  FOR  T-37  A/C 


SCHEDULE  3 


10 
20 
30 

40 
50 
60 
70 
80  * 

90  # 

100  * 
110  EXP 
0,0/. 1, 


*****#*#***************#****#********************************-********« 
4  4 

*  AIR-TRAFFIC  CONTROL  SIMULATION  * 

*  * 
****«*»***t***t***t**«»*»**«**#»***W**f***t*«********«»-*«*tt* 
* 

EXPONENTIAL  DISTRIBUTION 

FUNCTION  RN1 ,C24 

104/. 2,. 222/. 3,. 355/. 4,. 50?/. 5A. 69/, 64. 915/. 7, 1.2/. 75, 1.38 


W  •  1  ,  .  1  V*t/  .  Z  ,  .  ZZZ/  .  O  ,  .  OJJ/  .1..  •  U,  .  07/  .0..7IJ/  .  I  ,  1  •  z/  .  r  .  i  .  ..'e 

8. 1.6/. 84.1. 83/. 88, 2. 12/. 9,2. 3/. 92,2. 52/. 94, 2. 81/. 95,2. 99/. 96, 3. 2 
97, 3. 5/. 98,3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 


120  * 

130  * 

140  MEAN  FUNCTION  C1,DB 
3600 , 380 . 1 1 /5400 ,221.60/13200.488. 62 
1 5000 , 242 . 63/23400 ,557.83/ 25200 ,124.58 
31800, 494 . 60/36600 ,221.30 
150  * 


160  * 

290 
300 
310 
320 
330 
340  * 

350 
360 
370 
380 
390  * 

400  * 

410  * 

420 
430 
435 
450 
460 
470  * 

510  ACFTT2 

520  DOWN 

530 

540 

550 

560 

570  LINUP 
580 
590 
600 
610 
615 
620  # 

630 

640  AREA 
650  * 

660  * 

670  * 

680 

690 

700  WAIT 


GENERATE 

SEIZE 

ADVANCE 

RELEASE 

TERMINATE 

GENERATE 

ASSIGN 

MARK 

TRANSFER 


9000,FN*EXP 

DANGER 

600,180 

DANGER 


9100,FN$EXP 

4’2 
, TOEP 1 


; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

; A/C  FROM  OTHER  BASES 
•ASSIGN  TO  A  PARAMETER 
{MARK  THE  TIME 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


GENERATE 
TEST  L 
TRANSFER 
ASSIGN 
TRANSFER 

ASSIGN 

QUEUE 

GATE  FV 

GATE  SE 

SEIZE 

DEPART 

ADVANCE 

RELEASE 

MARK 

MARK 

ADVANCE 

TRANSFER 

SELECT  E 
TEST  E 


FN*M£ AN , FNf  £ X P 
N4GNRT ,  CP  I  ST  A 
.  49 , ,  ACFTT2 


l&OWN 


1,0 

CNTR 

DANGE 

CAP6 

CNTR 

CNTR 

70,20 

CNTR 

6 

7 

100,10 

.25,, OTHER 


2,101.108,0^.01 
B$*FlftST, i ,SEST 


QUEUP 


EAST  AREAS 


SAVEVALUE 

SAVEVALUE 

QUEUE 


13+,P1 
14+,  1 
P2 


; CREATE  AIRCRAFT 

•THE  SIMULATION  STOPS  AFTER 

; ASSIGN  TO  A  PARAMETER 

s ASSIGN  TO  A  PARAMETER 
{WAIT  OUT  OF  THE  RUNWAY 
•,DONT  MOVE  IF  EMERGENCY  HOLDS 
{CHECK  FOR  A\C  ON  BASE  LEG 
{CAPTURE  THE  CONTROLLER 
■GOING  FOR  LINE  UP 
•LINE-UP  CHECK 
•TAKE-OFF 

START  FLIGHT  TIMEFOR  T-2  A/C 
•START  FLIGHT  TIME  FOR  T-37  A/C 
{AFTER  TAKE  OFF  TO  DEPART  POINT 
{FORMATIONS,  INSTR.FL. ,  CPM 

{FIND  EMPTY  AREA  IF  EXISTS  ONE 
{fcHOOSE  EAST  OR  NEST  AREA 


{RECORDS  T-37 'S 
■RECORDS  TOTAL  #  OF  A/C 
{WAIT  IN  THE  QUEUE 
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APPENDIX  C 


This  appendix  contains  the  historical  data  for  a  three  day 
period.  For  each  aircraft  the  historical  data  consisted  of  the 
interarrival  times  for  take  off,  the  time  spent  waiting  in  the 
air  traffic  controller  queue,  the  time  spent  waiting  to  enter 
area  3  and  the  total  flight  time  of  T-37  aircraft.  All  times  are 
stated  in  seconds. 


Interarrival  times  for  take  off 


INARRDAY1  ! Interarrival  Tiles  for  tike  off  Diy  1) 

557  218  254  399  963  46  16  91  211  490  484  622  43  11  7  541  41  171  160  177  191  435 
472  62  27  646  305  1410  154  490  196  479  358  1637  207  84  568  548  1211  511 
31  225  65  565  41  149  778  279  157  362  575  1147  305  288  475  383  635  118  644 
274  1623  117  173  6  31  58  24  57  183  775  92  7  142  260  7  461  77  260  587  77 
39  378  225  235  766  70  183  200  79  322  22  59  118  298  70  3  280  146 


INARRDAY2  ( Interarri val  Ti»es  for  tike  off  Diy  2) 

291  228  259  597  89  70  351  1415  147  179  1006  388  167  43  478  22  135  18  21  357  257 
606  207  413  17  260  95s?  414  1B2  1424  486  33  213  304  1168  362  261  177  211 
91  53  310  14  128  267  304  90  37  274  11  804  959  580  425  1491  280  887  6  235 
947  114  155  141  157  46  63  152  155  151  113  128  37  1101  1142  919  242  132 
1460  216  397  61  1029  393  26  221  80  52  534  85  99  261  382  60  112 


INARRDAY3  (Interirrml  Tiies  for  tike  off  Diy  3) 

438  170  1353  577  124  144  88  338  260  636  330  27  547  15  427  52  127  42  267  23  11 
94  323  HI  596  115  798  421  148  171  657  214  20  371  106  1356  7B  120  22  51 
114  64  922  1135  481  1832  665  296  45  34  519  230  204  14  102  708  1238  1134 
1027  570  564  916  221  49  333  7  168  87  91  2  504  252  35  62  39  55  144  165  67 
107  2  116  146  357  270  980  1422  869  334  1097  873  277  447  312  621  274  30  21 
492  71  452  530  741  143  676 
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Time  spent  waiting  In  the  air  traffic  controller  queue 


Waiting  Controller  queue  Area  3 

Time  dayl  day2  day3  dayl  day2  day3 


0 

42 

45 

49 

14 

12 

13 

60 

42 

38 

31 

0 

0 

1 

120 

12 

20 

11 

1 

1 

0 

180 

3 

5 

4 

0 

1 

0 

240 

2 

2 

0 

0 

0 

1 

300 

0 

0 

0 

0 

1 

0 

420 

0 

2 

0 

0 

0 

0 

480 

0 

0 

0 

0 

0 

1 

540 

0 

0 

1 

0 

0 

0 

660 

0 

0 

0 

1 

0 

0 

720 

0 

0 

0 

1 

0 

0 

840 

0 

0 

0 

0 

0 

1 

960 

0 

0 

0 

1 

0 

0 

1020 

0 

0 

0 

0 

1 

0 

Total  flight  time  aircraft  T-37 


Time 

dayl 

day  2 

day3 

3960 

6 

5 

8 

4080 

23 

18 

13 

4200 

12 

23 

17 

4320 

10 

9 

8 

4440 

4 

3 

2 

4560 

0 

2 

1 

4800 

0 

1 

0 
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APPENDIX  D 


This  appendix  contains  the  summary  statistics  and  the 
hypothesis  test  results  for  the  eight  time  periods  for  which 
the  inter arrival  times  for  take  off  were  homogeneous.  These 
times  periods  were: 


1) 

07:30  - 

08:40 

5) 

11:40  - 

14:00 

2) 

08:40  - 

09:10 

6) 

14:00  - 

14:30 

3) 

09:10  - 

11:10 

7) 

14:30  - 

16:20 

4) 

11:10  - 

11:40 

8) 

16:20  - 

17:30 

The  Tables  D.l  and  D.2  contain  the  summary  statistics 
for  each  time  period  and  the  figures  D.l  through  D.8  contain 
the  distribution  fitting  and  the  Kolmogorov- Smirnov  test 


results. 


3  iO 


Table  D.l  Summary  Statistics 


Va  i  i  a  Lie: 

INTERVAL 1 

I NTERVAL2 

Sample  ai:e 

35 

45 

Average 

380.114 

221.689 

Median 

260 

160 

Mode 

259 

11 

G  t.-  o  fii  c?  t  r"  .i  c  m  t?  «  fi 

24B.624 

116.221 

Var lanes 

120438 

44135.7 

G  tandard  dev  i  a  tion 

347.041 

210.005 

S 1;  a  n  d  a  r  d  err  o  r 

58 . 6607 

31.3176 

Min  imum 

16 

7 

Ma:  l  imiiti 

1415 

798 

Range 

1399 

791 

Low  o  r  q  u  a  r  1.  i  1  e 

144 

42 

Upper  quantile 

547 

413 

lr.terqu.art)  It*  range 

403 

371 

L-wnusi 

1.65105 

0 . 928257 

Standardized  si ewness 

3 . 98766 

2.54214 

Kurtosi s 

2.56701 

-0.118849 

St  aridardi e d  1  u r t o s i s 

3 . 09996 

-0.162741 

Vai  table:  INTERVAL.?  I NTERVAL4 


-1  ll  i  p  3.  C.J  Z'j  1  II  0:' 

47 

27 

A  /  e?  r  3  q  c-2 

488.617 

242.63 

Madian 

305 

204 

Mode 

304 

14 

Geometric:  mean 

284.774 

144.641 

V  jfM^nc  p 

229164 

49612.2 

Standard  d  <  ?  v  i  a  t  i  o  n 

478.711 

222 . 738 

Standard  error 

69.8271 

42.866 

Min  imum 

20 

14 

Ma'  imum 

1832 

778 

Range 

1812 

764 

L  owe r  quart l  1 e 

120 

45 

Upper  quart i le 

646 

304 

Interquartile  range 

526 

259 

FT  ewness 

1.28883 

1 . 1451 

Standardized  skewness 

3 . 60734 

2.42911 

t urtosis 

0.696499 

0.434&59 

Standardized  kurtosis 

0.974683 

0.461025 
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Table  D.2  Summary  Statistics 


Viir  iabl  e : 

INTERVALS 

INTERVALS 

Sample  size 

3  6 

39 

Average 

557.833 

124.897 

Median 

450 

107 

Mode 

383 

Geometric  mean 

322 . 289 

69.1587 

Variance 

190050 

20045.8 

Standard  deviation 

435.947 

141.583 

Standard  error 

72.6579 

22.6714 

Minimum 

6 

Max  imu  m 

1623 

775 

Ra»nc|tt 

1617 

773 

Lower  quartile 

228 

39 

Upper  quartile 

901.5 

155 

Interquartile  range 

673 . 5 

116 

Skewness 

0.705837 

3.12258 

S tandardi zed  skewness 

1.72894 

7.96105 

Kur tosis 

-0 . 301 500 

12.1983 

S t  a n  ( J  a r  d i zed  k  u  r  t  a s 1 s 

-0 . 36927 

15.5498 

V  a  r  l  a  b  .1  e : 

INTERVAL? 

INTERVALS 

Sample  size 

38 

30 

Average 

494 . 605 

221.3 

Median 

345.5 

130.5 

Mode 

77 

112 

Geometric  mean 

327.136 

129.011 

Variance 

165600 

42845 

Standard  deviation 

406.94 

206 . 99 

Standard  error 

66.0143 

37.7911 

M j n  r mum 

26 

pr 

Maximum 

1460 

741 

Range 

1434 

738 

Lower  quartile 

216 

70 

Upper  quartile 

869 

Interquartile  range 

653 

Skewness 

0.90161 

1 . 10583 

Standardized  skewness 

2.269 

2.4727 

Kur tosis 

-0.289484 

0.246429 

Standardized  kurtosis 

-0.364259 

0.275516 

Figure  D.l  Distribution  Fitting  Interval  1 
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Estimated  KOLMOGOROV  statistic  DF’LUS  =  O. 
Estimated  KOLMOGOROV  statistic  DMINUS  =  O 
Estimated  overall  statistic  DN  =  0.134797 
Approximate  significance  level  =  0.3G6891 


134797 

.111456 


Figure  D.2  Distribution  Fitting  Interval  2 


Estimated  KOLMOGOROV  statistic  DF'LUS  = 
Estimated  KOLMOGOROV  statistic  DMINUS  = 
Estimated  overall  statistic  DN  =  0.0722 
Approximate  significance  level  =  1 


0 . 0708765 
0.0722221 


Figure  D.3  Distribution  Fitting  Interval  3 


Frequency  Histogram 


INTOL4 


Estimated  KOLMOGOROV  statistic 
Estimated  KOLMOGOROV  statistic 
Estimated  overall  statistic  DN 
Approximate  significance  level 


DPLUS  =  0.0935019 
DMINUS  =  0.0871474 
=  0.0935019 
=  1 


Figure  D.4  Distribution  Fitting  Interval  4 
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IHBML5 


Estimated  KOLMOGOROV  statistic  DPLUS  =  0.06380 
Estimated  KOLMOGOROV  statistic  DMINUS  =  0.1103 
Estimated  overall  statistic  DN  =  0.110322 
Approximate  significance  level  =  0.999951 


Figure  D.5  Distribution  Fitting  Interval  5 


r  j  xj 


Frequency  Histogram 


0  200  400  600  800 

mm 


Estimated  KOLMOGOROV  statistic 
Estimated  KOLMOGOROV  statistic 
Estimated  overall  statistic  DN 
Approximate  significance  level 


DPLUS  =  0 . 128466 
DM I NOS  =  0 . 0882622 
=  0.128466 
=  0.540454 


Figure  D.6  Distribution  Fitting  Interval  6 


Frequency  Histogram 


mm 


Estimated  KOLMOGOROV  statistic 
Estimated  KOLMOGOROV  statistic 
Estimated  overall  statistic  DN 
Approximate  siqnificance  level 


DPLUS  =  0.0867169 
DM I MUG  -  0.0760762 
=  0.0867169 
=  1 


Figure  D.8  Distribution  Fitting  Interval  8 
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APPENDIX  E 


This  appendix  contains  the  GPSS  program  outputs  for 
each  take  off  schedule. 

SCHEDULE  I 


FACILITY 

ENTRIES 

!  UTIL. 

AVE.  TIME 

AVAILABLE 

OWNER 

PEND 

INTER 

RETRY 

r-_ftY 

101 

14 

0.  877 

2130. 86 

1 

0 

0 

0 

0 

0 

102 

14 

O.E50 

2113.79 

1 

0 

0 

0 

0 

0 

103 

i  n 

4  4_ 

0.747 

216S.25 

1 

0 

0 

0 

0 

0 

104 

1  -> 

4  4. 

0.831 

2412.67 

1 

0 

0 

0 

c 

105 

11 

0.694 

2195.82 

1 

o 

0 

0 

0 

'*} 

106 

i  n 

4  *- 

0.80B 

2345.33 

1 

0 

0 

0 

0 

0 

107 

12 

0.799 

2317.33 

1 

0 

0 

0 

0 

0 

108 

9 

0.594 

2300.00 

1 

0 

0 

0 

0 

0 

DANGE 

o 

0.041 

726.00 

1 

0 

0 

0 

0 

0 

CNTR 

130 

0.266 

71.36 

1 

0 

0 

0 

0 

0 

DU^Yl 

116 

0.011 

3.35 

1 

0 

0 

c 

0 

0 

DUMY2 

160 

0.000 

0. 10 

1 

0 

0 

0 

0 

6 

DUMY3 

160 

0.000 

0.00 

1 

0 

0 

0 

0 

0 

DUMY4 

624 

0.047 

2.67 

1 

0 

0 

0 

0 

0 

QUEUE 

MAX 

CCNT. 

ENTRIES 

ENTRIES (0) 

AVE.C0NT. 

AVE. TIME 

AVE. (-0) 

RETRY 

101 

1 

0 

14 

11 

0.10 

255 . 07 

1190.33 

0 

102 

1 

0 

14 

13 

0.  03 

71.43 

1000.00 

fi 

103 

1 

0 

12 

11 

0.01 

29.92 

35 - . 00 

0 

104 

1 

0 

4  o 

4  42 

7 

0.04 

10- . 50 

262.80 

0 

1 05 

1 

0 

11 

9 

0.02 

53.00 

291.50 

0 

106 

1 

0 

i  n 

4  4- 

7 

0.01 

41.92 

100.60 

0 

107 

1 

0 

12 

7 

C .  05 

*  trcr  nc 

*  wU« 

372.60 

0 

108 

1 

0 

9 

5 

0.  04 

157.44 

354.25 

o 

CNTR 

"> 

4- 

0 

130 

-*v.‘ 

0.  09 

24.61 

32.93 

0 

CAP1 

*■> 

4. 

0 

116 

109 

0.01 

3.82 

63.29 

0 

CAD2 

1 

0 

160 

159 

0.00 

1. 10 

16.00 

0 

CAP3 

1 

0 

160 

160 

0.00 

0.00 

0.00 

0 

CAP4 

7 

0 

624 

557 

0.06 

3.11 

28.94 

0 

STORAGE 

CAP. 

REMAIN. 

MIN. 

MAX. 

ENTRIES 

AVL. 

AVE. C. 

UTIL. 

RETRY 

TELAY 

CAP1 

2 

2 

0 

2 

116 

1 

0.40 

0.201 

0 

0 

CAP2 

2 

o 

4. 

0 

2 

160 

1 

0.32 

0.161 

0 

0 

CAP3 

2 

2 

0 

*■> 

4- 

160 

1 

0.30 

0.150 

0 

0 

CAP4 

3 

3 

0 

3 

624 

1 

1.28 

0.425 

0 

0 

CAP5 

3 

3 

0 

3 

624 

1 

0.81 

0.269 

0 

0 

CAP6 

2 

2 

0 

2 

648 

1 

0.84 

0.418 

0 

0 

73 


A?!_E 

MEAN 

ST-.  DEO. 

RETRY  RAN3E 

FREQUENT' 

CNTR 

i.  ■» 

22.11 

r-, 

v_ 

- 

0 

77 

0  - 

60 

e4' 

60  - 

120 

1C 

CAP1 

3.82 

17.01 

0 

- 

0 

109 

30  - 

60 

4 

60  - 

90 

1 

90  - 

120 

1 

120  - 

150 

1 

CAP2 

0. 10 

1.26 

0 

— 

0 

159 

0  - 

30 

1 

CAP  3 

0.00 

0.00 

0 

— 

0 

160 

CAP4 

3.11 

9.92 

0 

— 

0 

557 

0  - 

30 

38 

30  - 

60 

29 

T0TAL37 

4159.52 

164.66 

0 

3840  - 

4080 

26 

4080  - 

4320 

31 

4320  - 

4560 

6 

4560  - 

4800 

T 

T0TA..T2 

4750.58 

138.22 

0 

4320  - 

4560 

4 

4560  - 

4800 

36 

4800  - 

5040 

21 

5040  - 

5280 

3 

TABLE 

MEAN 

STD. DEV. 

RETRY 

RANGE 

FREQUENCY 

AREA1 

255.07 

510.60 

0 

0 

11 

960 

- 

1080 

1 

1200 

- 

1320 

1 

1320 

— 

1440 

1 

AREA2 

71.43 

267.26 

0 

0 

13 

960 

- 

1080 

1 

AREA3 

29.92 

103.63 

0 

0 

11 

240 

- 

360 

1 

AREA4 

109.50 

181.31 

0 

0 

7 

0 

— 

120 

2 

240 

- 

360 

1 

360 

- 

480 

2 

AREA5 

53.00 

138.98 

0 

0 

9 

120 

- 

240 

1 

360 

— 

480 

1 

74 

Eur. 


25.38 

82.31 


93.97 

97.4! 

98.28 

99.14 

100.00 


99.77 

100.00 


100.00 


88.28 
95.35 
100. 00 


39.3? 

86.36 

95.45 

100.00 


6.25 

62.50 

95.31 

100. 00 


CUM.  •/. 


78.57 

85.71 

92.86 

100.00 


92.86 

100.00 


91.67 

100.00 


58.33 
75.00 

83.33 
100.00 


81.82 

90.91 

100. Cl 


I 


TAB_E 

AK'EAf: 

AREA? 


AREA8 


MEAN 
a;  .  92 


Sir. DEV.  RETRY  RANGE 


c-pr«-  r> 


63.52 


~7~? 


157.44 


283.67 


0 

120 


0 

240 

360 

600 


0 

120 

840 


0 

120 

240 


0 

120 

360 

480 

720 


240 

960 


7 

1 

2 
T 
i 


i 

n 


r;  —  — 


BZ.ZZ 


58.33 

66.67 

83.33 

91.67 

100.00 


55.56 
66 .  6 ; 
88.89 

100. oc 


1A 


FACILITY 

ENTRIES 

JTI ; 

AVE.  TINE 

AVAILABLE 

OWNER 

-'EM 

INTER 

SE~En 

BE- 

101 

4  ~T 

i  O 

0.  775 

2066.31 

1 

0 

0 

0 

0 

0 

102 

13 

0.792 

2111.69 

1 

0 

0 

0 

0 

6 

103 

11 

0.740 

2330. 18 

1 

0 

0 

0 

0 

V.‘ 

104 

11 

0.735 

2316.09 

1 

0 

0 

0 

0 

0 

105 

9 

0.600 

2311.22 

1 

0 

0 

0 

0 

0 

106 

9 

0.607 

2338.22 

1 

0 

0 

0 

0 

0 

107 

6 

0.417 

2410.33 

1 

0 

0 

0 

0 

c 

108 

6 

0.395 

2281.83 

1 

0 

o 

0 

0 

0 

DANSE 

A 

0.022 

772.00 

1 

A 

0 

0 

0 

0 

0 

CNTR 

106 

0.208 

67.97 

1 

0 

0 

0 

0 

6 

DUMYI 

93 

o.oi: 

4.33 

1 

0 

0 

6 

0 

0 

DUMY2 

123 

0.000 

0.11 

1 

0 

0 

0 

0 

0 

DUMY3 

123 

0.000 

0.  12 

1 

0 

0 

0 

0 

0 

DUMY4 

522 

0.019 

1.28 

1 

0 

0 

0 

0 

f) 

QUEUE 

MAX 

CONT. 

ENTRIES 

ENTRIES (0)  AVE. CONT. 

AVE. TIME 

AVE. (-0) 

RETRY 

101 

1 

0 

13 

13 

0.00 

0.00 

0.  on 

0 

102 

1 

0 

13 

13 

0.00 

0.00 

0.00 

0 

103 

1 

0 

11 

10 

0.00 

3.45 

38.  00 

fi 

10A 

1 

0 

11 

10 

0.00 

0. 18 

2.0C 

0 

105 

1 

0 

9 

9 

0.00 

0.00 

0.00 

0 

106 

1 

0 

9 

B 

0.01 

26.56 

239.00 

0 

107 

4 

i 

0 

6 

6 

0.00 

0.00 

0.00 

0 

ICE 

l 

0 

t 

C T, 

0.00 

13.00 

7S .  or 

o 

CNTR 

1 

A 

0 

106 

78 

0.01 

2.72 

10.2? 

c 

CAP1 

0 

93 

87 

0.01 

4.71 

73.00 

0 

CAP  2 

0 

123 

122 

0.00 

0.11 

14.  00 

o 

CAP  3 

1 

(': 

4  n-f 

A  a---' 

122 

0.00 

0.  12 

15.00 

0 

CAP4 

2 

0 

522 

494 

0.02 

1.33 

24.82 

0 

STORAGE 

CAP. 

REMAIN* 

MIN. 

MAX. 

ENTRIES 

AVL. 

AVE. C. 

UTIL. 

RETR", 

DELAY 

CAF'l 

O 

2 

0 

A, 

93 

1 

0.32 

0.  162 

0 

0 

CAP  2 

A, 

2 

0 

n 

Am 

123 

1 

0.25 

0. 123 

0 

0 

CAP3 

O 

n 

0 

n 

Am 

123 

1 

0.23 

0.  lit 

0 

'“l 

CAP  4 

7 

0 

*7 

522 

1 

1.06 

0.355 

0 

0 

CAP5 

3 

3 

0 

3 

522 

1 

0.67 

0.225 

0 

o 

CAP6 

n 

n 

A. 

0 

n 

Am 

572 

1 

0.74 

0.372 

6 

r 

TABLE 

MEAN 

STD. DEV. 

RETRY  RANGE 

FREQUENCY 

CUM.  7. 

CNTR 

2.72 

5.60 

0 

— 

0 

78 

73.58 

0  - 

60 

28 

100.00 

CAP1 

4.71 

18.51 

0 

— 

0 

87 

93.55 

30  - 

60 

2 

95.70 

60  - 

90 

oo.  t?2 

90  - 

120 

\ 

i no \ 

76 


i 


▼cpee 

Mrjv 

C tt\  nr\  ■ 

RETC'Y 

RANGE 

FRE'2_'EN3V 

p-  ^ 

c . : : 

1 . 

0 

- 

30 

*~1 

100. 00 

CAP™ 

0. 12 

<  TC 

X  «  -‘O 

0 

/■ 

1  O'-! 

C3i  ■  C 

0 

- 

30 

1 

100. O' 

CAP4 

1.33 

6.27 

0 

_ 

0 

494 

94. 6  A 

0 

— 

30 

19 

98.28 

30 

60 

9 

100.00 

T0TAL37 

4128.00 

136.31 

0 

3840 

- 

4080 

25 

48.08 

4080 

- 

4320 

on 

4.X. 

90.33 

4320 

— 

4560 

5 

1 00 . 00 

TDTALT2 

4703 . 52 

126.07 

0 

4320 

4560 

6 

11.11 

4560 

- 

4800 

37 

79.63 

4800 

- 

5040 

9 

96.30 

5040 

- 

5280 

100.00 

TAE^E 

MEAN 

STD. DEV. 

RETRY  RANGE 

FREQUENCE 

cut.  ; 

AREA1 

0.00 

0.00 

0 

0 

13 

100.00 

AREA2 

0.00 

0 . 00 

0 

- 

0 

13 

100.00 

AREA3 

3.45 

11.46 

0 

0 

10 

90.91 

0  — 

120 

1 

1 0  v1 .  0  -.: 

AREA4 

0.  18 

0.60 

0 

- 

0 

10 

9C.91 

o  - 

120 

1 

100.0 

A°EA3 

C  .  00 

0.00 

0 

0 

9 

100.00 

AREA6 

26.56 

79.6™ 

0 

0 

8 

88.  8° 

120  - 

240 

X 

100. 00 

AREA7 

0.00 

0.00 

0 

- 

0 

6 

100.00 

AREA8 

13.00 

<r 

00 

*— « 

K> 

0 

- 

0 

5 

83.33 

0  - 

120 

1 

100.00 
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SCHEDULE  IE 


FAC  11- 1 TY 

ENTRIES  UTIL. 

AVE.  TIME 

AVAILABLE 

OWNER 

PEND 

INTER 

RETRY 

DELA 

101 

12 

0.774 

2250. 17 

1 

0 

0 

0 

0 

0 

102 

10 

0.626 

2182.40 

1 

0 

0 

0 

6 

103 

9 

0.528 

2045.89 

1 

0 

0 

0 

0 

6 

104 

9 

0.549 

2127.67 

1 

0 

0 

0 

0 

f. 

105 

9 

0.578 

2238.22 

1 

0 

0 

0 

0 

0 

106 

c 

C.  619 

2399.67 

1 

c 

0 

0 

0 

107 

9 

0.637 

2466.11 

1 

0 

0 

0 

0 

c 

ioe 

6 

C.  541 

2357.25 

4 

1 

0 

0 

0 

0 

DANSE 

6 

0.097 

568.50 

1 

0 

0 

0 

0 

0 

CNTR 

89 

0. 182 

71.39 

1 

0 

0 

0 

0 

0 

DUMY1 

63 

0.014 

8.25 

1 

0 

0 

0 

o 

0 

DUMY2 

9B 

0.002 

0.92 

4 

A 

0 

0 

0 

/v 

DUMY3 

98 

0.000 

0.00 

1 

0 

0 

0 

6 

0 

DUMY4 

423 

0.  009 

0.81 

4 

A 

0 

0 

0 

0 

c 

QUEUE 

MAX 

CONI. 

ENTRIES 

ENTRIES (0) 

AVE.CONT. 

AVE. TIME 

AVE. (-0; 

RETRY 

101 

1 

0 

12 

8 

0.09 

251.83 

755.50 

0 

102 

1 

0 

10 

9 

0.03 

108.80 

1088.00 

0 

103 

1 

0 

9 

9 

0.00 

o .  or- 

c.oc 

/* 

104 

1 

0 

9 

9 

0.00 

0.00 

0.  00 

0 

105 

1 

0 

9 

9 

0.00 

0.00 

0.00 

r\ 

K' 

106 

1 

0 

9 

9 

0.00 

0.00 

0.00 

0 

107 

1 

0 

9 

9 

0.00 

0.00 

0.00 

0 

108 

1 

0 

8 

8 

0.00 

0.00 

o.oc 

o 

CNTR 

'n 

0 

89 

24 

0.07 

25. 52 

3  C .  9  C 

V 

CAP1 

n 

0 

63 

54 

0.02 

9.57 

67 . 00 

o 

CAP2 

i 

6 

98 

96 

0.00 

0.92 

45 . 00 

c 

CAP3 

i 

0 

98 

98 

0.00 

0.00 

0 . 00 

CAP4 

n 

o 

423 

410 

0.01 

0.85 

27.62 

o 

STORAGE 


CAP1 

CAP2 

CAP3 

CAP4 

CAP5 

CAP6 


CAP.  REMAIN'. 


2  2 

2  2 

2  2 

3  3 

3  3 

2  2 


MIN. 

MAX. 

ENTRIES 

AVL. 

AVE. C. 

UTIL. 

RETRY  DE-£ 

0 

2 

63 

1 

0.22 

0.110 

0  0 

0 

2 

98 

1 

0.20 

0.099 

0  0 

0 

2 

98 

1 

0.18 

0.092 

0  0 

0 

3 

423 

1 

0.86 

0.286 

0  0 

0 

3 

423 

1 

0.55 

0.184 

0  0 

0 

2 

481 

1 

0.62 

0.308 

0  c 

78 


~A3lE 

CK'TC 

cap: 

CAP2 

CAP3 

CAPA 

T4B,_E 

T07AJ37 

T0TALT2 


TABLE 

AREA1 

AREA2 

AREA3 

AREA4 

AREAS 

AREA6 

AREA? 

AREA3 


ST3.DEV.  RETRv  t?£\i£r 


o 


0.92 


0.85 


4680.07 


r  - 

60  - 


30  - 
60  - 
90  - 

7.87  0 

0  - 
60  - 

C .  v  0 


60 


0 

60 

90 

120 


0 

30 

90 


0 


5.24  0 


0  - 
30  - 


Rc.TRV  RAN3E 


0 

30 

60 


1 6B . 33  0 

3840 

4080 

4320 

4560 

103.17  0 

4320 

4560 

4800 


4080 

4320 

4560 

4800 


4560 
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